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ABSTRACT 


Energy  transfer  in  chemically  reacting  boundary  layer  flows  is  discussed 
from  the  point  of  view  of  the  investigator,  who  is  seeking  to  extend  existing 
correlation  formulae  to  cases  in  which  thermochemical  effects  influence  heat 
transfer  rates.  Emphasis  is  placed  on  the  prediction  of  convective  heat  fluxes 
in  high  performance  rocket  motors;  however,  examples  are  also  taken  from  the 
field  of  hypersonic  gas  dynamics. 

The  following  topics  are  cons'  '-r 

1.  the  appropriate  driving  force  for  heat  transfer  with 
chemical  reaction 

2.  effects  of  the  enhanced  efficiency  of  energy  transport 
by  diffusion  as  compared  to  ordinary  conduction 

3.  calculation  of  the  turbulent  film  conductance  in  axi- 
symmetric  nozzles 

4.  thermodynamic  calculation  of  enthalpy/mixture-ratio 
charts  for  combustion  gas  mixtures 

5.  effects  of  chemical  non-equilibrium  in  the  gas  phase 
(during  the  expansion  process  as  well  as  within  the 
boundary  layer) 

6.  effects  of  surface  catalyzed  exothermic  recombination 
reactions 

7.  estimation  of  transport  properties  in  partially  disso¬ 
ciated  gas  mixtures  with  emphasis  on  the  binary  diffu¬ 
sion  coefficients  pertaining  to  molecular  fragments. 

Several  detailed  calculations  are  included  for  the  case  of  pure  dissociating 
hydrogen  and  for  the  products  of  combustion  of  hydrogen  and  oxygen.  Areas  in 
need  of  additional  investigation  are  pointed  out,  and  extensive  references  are 
made  to  recent  work  which  is  felt  to  have  a  bearing  on  the  topics  discussed. 
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exponent  on  the  Prandtl  number  in  a  correlation  equation  of  the 
Nusselt  form 

catalytic  parameter  [eq.(7)] 
mass  fraction 

heat  capacity  (per  gram)  at  constant  pressure 

heat  capacity  (per  gram  mole)  at  constant  pressure 

heat  capacity  (per  gram  mole)  at  constant  volume 

local  nozzle  diameter 

binary  (Pick)  diffusion  coefficient 

mixture  ratio  parameter  [eq.(S6)] 

mass  velocity  (mass  flow  per  unit  area) 

static  enthalpy  of  mixture  (including  chemical  contributions) 
molar  enthalpy 

energy  flux  vector  relative  to  mass  averaged  velocity 
interfacial  reaction  rate  constant  (velocity)  [eg. (4)] 

Boltzmann  constant 
characteristic  physical  length 
chemically  "frozen"  Ls>’"'is -Semenov  number 

-  D/[X^/(pCp^ ^)] 

atomic  or  molecular  mass 
molecular  weight 

non-dimensional  (heat  or  mass)  transfer  coefficient  (Nusselt  number) 
total  number  of  components  in  mixture 
tota 1  pressure 
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Pr*  Pvandtl  number  for  heat  conduction  =  (|ji/p)/[X/ (pc  )] 

A  P 

Prjj  Prandtl  number  for  diffusion  =  (p,/p)/D 

(often  called  the  Schmidt  number  In  the  western  literature) 

4  heat  flux  (per  unit  area)  at  surface 

Q  heat  of  recombination  -  hi  -  hj  (per  gram) 

recovery  factor  tor  directed  kinetic  energy 

Tp  recovery  factor  for  chemical  energy 

r  mixture-ratio  parameter  (oxyqen-to-fuel  mass  ratio) 

R  universal  gas  constant  [  -  1.9872  ca 1  (qm-mole)  (  K)  j 

Re  Reynolds  number 

s  Reynold's  analogy  factor  (ratio  of  l/2  skin  friction  coefficient  to 

Stanton  number 

St  Stanton  number  (non-dimensional  heat  or  mass  transfer  coefficient) 

T  absolute  temperature  (usually  '^K) 

Uiv  X  and  y  component,  respectively,  of  the  local  gas  velocity 

x,y  physical  coordinates  parallel  and  perpendicular  (respectively)  to 

the  gas-solid  Interface!  x  Is  also  >ised  for  mole  fraction 

7.  stretched  streamwlse  distance  containing  aerodynamic  as  well  as 

chemical  kinetic  parameters 

p  Invlscld  velocity  gradient  at  the  nose  [eq.(3)] 

E  energy  parameter  appearing  In  Lennard-Jones  Interaction  potential 

function 

Y  recombination  probability  or  ratio  of  specific  heats 

r)  normalized  diffusionai  reduction  in  reactant  concentration  at  inter¬ 

face  [eo. (lO)] 

9  function  defined  in  text  [eg. (77)] 

\  thermal  conductivity  of  mixture 


X 


|ji  absolute  viscosity  of  mixture 

p  absolute  density  of  mixture 

<3  molecular  size  parameter  (in  angstrom  units) 

t:  characteristic  time 

<p  heat  flux  potential  [eq.(32)] 

<5  equivalence  ratio  I-  (Ha/0,)/(H,/0a)  .  .  , 

*  •  stoich 

function  defined  by  eq.(l03) 

function  of  chemical  kinetic  parameters  [see,  for  example,  eqs.(90) 

^  and  (97)] 

>1^1, 4^a  functions  defined  In  text  [see  eq.(50)] 

(d 

u  parameter  in  the  viscosity-temperature  law  p  ~  T 

qrad  gradient  operator 

A  change  in  (across  the  boundary  layer j  or  across  a  chemical  reaction, 

depending  on  uoiiLexl) 

d  ordinary  differential  operator 

fl  partial  differential  operator 

summation  over  i 

lim  limit  operation 


Subscripts 

1  pertaining  to  the  lighter  constituent  (atoms  in  the  case  of  a 
dissociated  diatomic  qas)  except  when  used  with  the  function  ijr 

2  pertaining  to  the  heavier  constituent  except  when  used  with  the 
function  \l/ 

A  pertaining  to  atoms  in  a  dissociated  diatomic  gas  or  one  of  the 

propellants  in  a  rocket  motor 

pertaining  to  the  difference  between  A  and  M  or  the  interaction 
between  A  and  M  depending  on  context 


y  1 


A-M  or 
AM 


avg 


average 


B  one  of  the  propellants  in  a  rocket  motor 

c  chamber  upstream  of  rocket  nozzle 

chem  chemical 

D  pertaining  to  diffusion 

e  at  outer  edge  of  boundary  layer 

eq  pertaining  to  thermochemical  equilibrium 

form  pertaining  to  the  "formation"  reaction 

f  chemically  frozen 

1  pertaining  to  species  i  or  pertaining  to  the  injector  station  in  a 

rocket  thrust  chamber 

j  pertaining  to  species  j 

kin  pertaining  to  ordered  kinetic  energy 

M  pertaining  to  molecules  in  a  dissociating  diatomic  gas 

m  pertaining  to  the  mixture 

r  recovery  (adiabatic  wall) 

stoich  stoichiometric 

T  thermal  (sensible)  contribution 

t  throat  section  of  nozzle 

w  at  the  wall  (surface  of  body) 

X  based  on  streamwise  distance  along  interface 

X  pertaining  to  thermal  conduction 

*  evaluated  at  the  reference  temperature  (or  reference  enthalpy) 

or  pertaining  to  some  critical  (singular)  value 
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I  lOTRODUCTORY  REMARKS 


A  goal  of  heat  and  mass  transport  theory  is  to  arrive  at  vyorking  for¬ 
mulae,  having  a  rational  basis,  which  can  be  used  to  predict  heat  transfer 
rates  to  solids  under  conditions  for  which  there  is  no  direct  experimental 
data  (i .0. ,  for  purposes  of  extrapolation  and  interpolation).  By  a  rational 
basis  we  imply  that  the  relations  are  consequences  of  a  realistic  physical 
model  which  draws  upon  the  interaction  of  more  or  less  elementary  processes, 
each  of  which  is  quantitatively  understood  as  a  result  of  Independent  investi¬ 
gation.  No  more  of  these  elementary  processes  are  to  be  invoked  than  are 
necessary  to  account  for  the  (necessarily)  limited  experimental  data  that  are 
available  at  any  given  time. 

Experience  in  seemingly  diverse  fields  has  taught  us  that  chemical 
change  can  cause  dramatic  effects  on  heat  transfer  rates  both  in  the  absence 
and  presence  of  convective  fluid  motion.  While  earlier  cases  can  be  cited, 
a  particularly  beautiful  example  is  provided  in  the  work  of  I.  Langmuir*^ 
(1912),  who  noticed  that  "at  extremely  high  temperatures  the  power  consump¬ 
tion  necessary  to  maintain  a  tungsten  wire  (Nernst  filament)  at  a  given  tem¬ 
perature  in  (gaseous)  hydrogen  increases  abnormally  rapidly  with  the  temper¬ 
ature".  What  followed  was  a  convincing  guantitatlve  demonstration  that  this 
"abnormal"  increase  was  due  to  the  "normal"  endothermic  fragmentation  of  the 
hydrogen  molecules  in  contact  with  the  filament,  followed  by  their  diffusion 
and  subsequent  exothermic  reassociation  (recombination)  in  the  cooler  regions 
of  the  gas  away  from  the  filament. 

The  present  work,  devoted  to  convective  heat  transfer  with  chemical 
reaction,  is  motivated  by  the  fact  that  propulsion  technology  has  attained 
a  state  of  development  which  obliges  us  to  consider  thermochemical  effects 
of  the  type  described  abovej  both  for  the  prediction  of  heat  transfer  rates 
to  external  surfaces  of  high  speed  veiiicles  and  internal  surfaces  (e, q. , 
the  prediction  of  roolinq  reauirements  for  high  performance  thrust 
chambers^  ^  3  J, 67, 90 
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In  attempting  to  understand  the  convec  ve  heat  exchange  phenomena 
occuring  when  solid  surfaces  interact  with  moving  fluids  with  which  they  are 
not  in  chemical  as  well  as  thermal  equlllbrlumi  it  is  to  be  expected  that  the 
number  of  independent  parameters  required  to  adequately  describe  the  ensuing 
processes  will  be  large.  For  the  first  time,  chemical  kinetic  properties, 
thermochemical  properties,  as  well  as  diffusive  properties  of  the  constituent 
fluids,  are  matters  of  great  consequence.  In  seeking  compact  genera lixatlons , 
we  are  demanding  more  and  more  since  correlation  formulae  must  predict  the 
unfamiliar  as  well  as  reduce  to  the  familiar.  As  a  result,  for  many  problems 
of  current  practical  interest,  it  is  not  a  priori  evident  that  a  reasonable 
degree  of  generality  can  be  achieved  without  excessive  sacrifice  of  accuracy. 

This  is  an  area  in  which  scientific  taste  and  temperament  seem  to  vary 
widely  and  one  which  will  profit  by  two  parallel  approaches.  On  the  one  hand, 
there  are  tiiosc  investiaa tors  who  will  point  to  the  conservation  equations  of 
aero-thermochemistry  (see,  for  example,  reference  59)  and  assure  the  designer 
that  any  particular  problem  he  wishes  solved  can  he  handled  by  machines  pro¬ 
vided  the  requisite  input  data  Is  available.  On  the  other  hand,  there  are 
those  investigators  for  whom  this  position  constitutes  a  sort  of  intellectual 
defeat,  the  argument  being  that  if  each  new  situation  must  be  handled  indepen¬ 
dently  there  is  no  limit  to  the  number  of  problems  which  will  have  to  be 
addressed  to  computing  machines,  nor  will  there  be  a  limit  to  the  space  re¬ 
quired  tn  present  the  results  to  the  scientific  community.  The  goal  of  this 
second  group  is  a  series  of  compact  generalisations  based  on  physical  and 
mathematical  "models"*  i .e. ,  working  corielation  formulae  which  will  display 
fiinctional  dependences  and  on  which  further  predictions  can  be  based,  even 
for  situations  which  are  apparently  dissimilar  to  those  for  which  the  gcneral- 
Irations  have  been  developed.  Only  theoreticians  in  the  first  group  and  ex¬ 
perimentalists  will  be  ablp  to  provide  the  yardsticks  necessary  to  judge  the 
relative  merits  (physical  reality)  of  each  approximation  or  model.  On  the 
other  hand,  theoreticians  in  the  second  group  can  provide  approximate  simili¬ 
tudes  which  are  of  conceptual  as  well  as  engineering  value.  Thus  the  interplay 
between  these  three  groups  will  determine  trie  rate  at  wnii,}i  oui  understanding 
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will  improve 


A  designer  consulting  the  literature  for  direction  rarely  has  the  good 
fortune  of  finding  experimental  data  on  the  very  system  in  question  over  the 
appropriate  range  of  variables,  nor  has  he  the  time  or  money  to  adopt  cut-and- 
try  methods.  As  a  result,  if  rapid  advances  are  to  be  made  there  is  no  alter¬ 
native  but  to  generalize  existing  information  to  the  greatest  extent  possible 
and  to  express  the  results  in  a  useable  form.  To  merely  learn  that  a  particular 
problem  is  in  principle  solved  Is  usually  sm=ill  consolation.  We  are  reminded 
of  the  words  of  H.  Poincare^ t  "It  is  far  better  to  predict  without  certainty, 
than  never  to  have  predicted  at  all". 

In  what  follows,  the  subject  of  energy  transfer  in  chemically  reacting 
boundary  layer  flows  is  discussed  from  the  point  of  view  of  the  investigator 
who  is  seeking  to  extend  existing  heat  transfer  correlation  formulae  to  cases 
in  which  thermochemical  effects  influence  heat  transfer  rates.  While  the 
present  discussion  will  primarily  bo  directed  at  the  prediction  of  rocket 
motor  heat  fluxes,  examples  taken  from  the  field  of  hypersonics  will  also  be 
included  when  these  are  felt  to  shed  additional  light  on  the  class  of  pheno¬ 
mena  being  discussed. 

The  cooling  prolilem  is  important  enough  to  be  the  limiting  factor  in 
the  design  of  many  compact,  hinh-performance  thrust  chambers.  With  the  use 
of  more  eneroetic  chemical  propellant  combinations  [  i  .e. ,  propellant  combin¬ 
ations  yielding  higher  values  of  the  characteristic  velocity  c"'-(T^A^)^  ]  the 
convective  heat  transfer  rates  everywhere  within  the  chamber  will  increase. 

An  important  question  is  the  extent  of  this  increase.  It  is  for  these  pro¬ 
pellant  combinations  that  the  thermochemical  and  diffusion  effects  to  be  dis¬ 
cussed  are  likely  to  be  most  noticeable,  since  a  substantial  fraction  of  the 
combustior  products  at  the  chamber  temperature  and  pressure  are  in  the  form  of 

^Science  and  Hypothesis.  Hypotheses  in  Physics,  Chapter  IX,  Part  IV,  p. 144, 
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light)  dissociated  gases. 


Convective  heat  transfer  in  rocket  motors  has  already  been  discussed 
from  many  points  of  view  (see,  for  example,  references  6,7,22,24.41,115,116, 
117,129)  but  few  treatments  have  Included  more  than  several  cursory  remarks 
on  the  subject  of  thermochemical  effects.  It  should  not  be  implieci  that  a 
sufficient  amount  of  data  and  theory  has  accumulated  in  the  interim  to  enable 
accurate  quantitative  predictions  in  this  difficult  area.  But,  on  the  other 
hand,  it  is  felt  that  the  field  has  progressed  on  some  fronts  beyond  the  point 
of  qualitative  speculation.  It  is  still  true  that  in  actual  liquid  propellant 
motors  the  effect  of  changing  the  injection  pattern  may  far  outweigh  several 
of  the  effects  to  be  described.  But,  this  fact  alone  should  not  be  allowed 
to  hinder  the  development  of  the  theory  of  convective  heat  transfer  with  chem¬ 
ical  reaction.  On  the  contrary,  the  more  accurate  our  Ideali^'ed  predictions 
become,  the  more  we  will  be  able  to  say  about  the  actual  effects  of  Injection 
pattern  and  other  contributions  to  the  energy  flux  (e.q, ,  radiation)  in  the 
future. 

The  present  work  represents  an  attempt  to  provide  some  of  the  answers 
to  the  questions  posed  above.  Among  the  topics  discussed  are:  the  driving 
force  for  heat  transfer  with  chemical  reactions  factors  accoutilliig  for  the 
enhanced  efficiency  of  energy  transport  by  diffusions  the  calculation  of  the 
turbulent  film  conductance  in  axi-symmetrlc  noxxless  the  calculation  of  thermo¬ 
dynamic  charts  for  propellant  gases.  In  pa i  ti c.ulfir .  hydrogen  and  oxygen  com¬ 
bustion  products;  the  effects  of  chemical  non-equilibrium  in  the  gas  phase 
during  the  expansion  process  and  within  the  boundary  layer,  as  well  as  the 
effects  of  the  surface  catalyzed  exothermic  atom  recombination;  and  lastly, 
the  estimation  of  transport  properti«'s.  Particular  attention  will  be  paid  to 
the  case  in  which  the  mass  dlffusivity  for  leactive  species  is  different  from 
the  thermal  diffuslvity  of  the  gas  mixtures  i .e. ,  the  case  of  Lewi s-Semenov 
number  different  fiom  unity. 


/. 


II  THE  DRIVING  FORCE  FOR  CONVECTIVE  HEAT  TRANSFER 
WITH  CHEMICAL  REACTION 


In  the  study  of  transfer  problems  (heat  or  mass)  the  conceptual  pattern’^* 
followed  may  be  broken  down  as  followsi 

1.  a  "driving  force"  Is  defined  which  is  regarded  as  the  cause 
of  the  transfer 

2.  a  "coefficient"  is  defined  as  the  rate  of  transfer  per  unit 
driving  force 

In  this  scheme  the  values  of  these  coefficients  as  well  as  the  nature  of  the 
driving  force  become,  in  themselves,  legitimate  objects  of  study.  It  will  be 
appreciated  that  an  arbitrariness  evists  since  any  "driving  force"  may  be 
selected  provided  the  corresponding  "coefficient"  is  experimentally,  or  theo¬ 
retically,  determined.  This  arbitrariness  is  reduced  by  adopting  the  principle 
that,  among  the  various  possible  effective  driving  for'-es,  one  should  select 
those  which  impart  to  the  resulting  coefficients  the  greatest  generality 
(applicability  over  the  widest  possible  range  of  experimental  conditions). 
Equivalently,  one  attempts  a  kind  of  "separation  of  variables"  such  that  the 
effects  of  changes  in  physical  parameters  are  preferably  confined  to  either 
the  coefficient  or  the  driving  force  but  not  reflected  in  both.  Ihls  principle, 
of  course,  is  not  peculiar  to  transfer  theory,  but  is  common  to  all  the  sciences. 

In  the  field  of  forced  convection  heat  transfer,  perhaps  the  best  known 
example  of  this  separation  of  variables,  is  offered  by  Newton's  "law"  of  cool¬ 
ing!  l.e.,  the  statement  that  the  heat  flux  ^  should  be  proportional  to  the 
temperature  difference  AT  between  the  fluid  and  the  surface.  While  it  is  not 
suggested  that  the  proportionality  "constant"  h  (the  heat  transfer  coefficient) 
is  truly  constant,  It  ^  Implied  that  h  has  no  dependence  on  AT  Itself.  That 
Is  1.0  say,  h  should  reflect  only  changes  in  fluid  dynamic  parameters. 

In  the  presence  of  chemical  energy  release,  it  will  be  seen  that 
Newton's  "law"  of  cooling,  as  such,  ceases  to  be  useful  in  identifying  the 
appropriate  driving  force  for  heat  transfer  and  heat  transfer  coefficient! 
i . e. ,  even  In  the  absence  of  viscous  heating  the  energy  flux  will  no  longer  be 
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simply  proportional  to  tha  temparature  difference  AT  between  the  fluid  In  the 
free  stream  and  the  Interface^,  This  Is  perhaps  easiest  to  visualize  for  the 
case  In  which  enerqy  Is  also  transported  to  the  wall  by  the  diffusion  of  atoms* 
which  react  exothermically  at  the  wall  Itself  but  not  within  the  gaseous  boun¬ 
dary  layer.  Here  the  total  energy  transport  per  unit  time  and  area  is  given 
by  the  sum  of  two  terms,  and  The  first  of  these  terms  represents  the 
conductive  (convective)  contribution  to  the  heat  flux  and  is  proportions  1  to 
the  temperature  gradient  (in  the  fluid)  established  normal  to  the  Interface. 

The  second  of  these  la  the  contribution  due  to  thermo-chemical  energy  transport 
through  the  boundary  layer  principally  by  concentration  (pick)  diffusion, 
followed  by  exothermic  chemical  reaction  at  the  Interface.  While  the  temper¬ 
ature  difference  AT  Is  approximately  the  driving  force  for  the  contribution 
the  atom  (reactant)  concentration  difference  Ac^  is  approximately  the 
driving  force  for  with  the  result  that  the  sum  Is  neither  propor¬ 

tional  to  AT  alone  nor  Ac^  alone  but  is  determined  by  some  combination  of  these 
pa rameters. 

Chemically  Frozen  Boundary  Lavers  with  Catalytic  Surface  Reaction 

If>  the  field  of  high  speed  flight,  this  particular  problem  has  been 
examined  for  the  case  of  laminar  stagnation  point  heat  transfer  to  blunt-nosed 
bodies  In  the  presence  of  catalytic  surface  reaction.**’**  Adopting  a  straight¬ 
forward  simill,tude  approach.  It  is  possible  to  derive  an  instructive  correl¬ 
ation  formula  for  the  heat  flux  which  can  then  be  used  as  the  starting  point 
for  a  more  general  discussion.  The  development  here  initially  parallels  that 
given  in  reference  (92). 


hlRt.orlr.al  Interest,  one  may  cite  the  papers  of  Rocard  and  Veron**'****’’' 
in  which  exothermic  chemical  change  Is  said  to  contribute  a  "convection  vlve" 
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We  write  the  contribution  to  ordinary  thermal  conduction  (convection) 
as; 

^  Stj^  G  Ah^  (1) 

where  Zih^  represents  the  difference  in  the  "frozen*'  (thermal)  specific  enthalpy** 
of  the  qas  at  the  outer  edge  of  the  boundary  layer  and  at  the  wall  Stj^  is  the 
non-dimensional  heat  transfer  coefficient  (Stanton  number),  and  G  is  the  mass 
velocity  p  u  .  The  energy  transfer  contribution  due  to  atom  diffusion  to 
the  wall  may  be  taken  as  the  product  of  the  heat  of  recombination  Q  (assumed 
constant)  and  the  rate  of  convective  diffusion; 

4^=0  Stp  G  Ac^  (2) 

The  sum  4  -  of  these  two  rates  will  be  the  observed  heat  transfer  rate. 

From  similitude  theory,  we  assume  that  tlie  Stanton  number  Stj^  for  mass  trans¬ 
port  by  convective  diffusion  is  obtainable  from  the  Stanton  number  Stj^  for 
heat  transport  by  making  the  replacement  Pr^^  Pr^  (Pr^  is  commonly  referred 
to  as  the  Schmidt  number  in  western  literature).  This  procedure  is  asympto¬ 
tically  exact  for  the  constant  property  case  when  the  free  stream  reactant 
concentration  (mass  fraction)  is  small  compared  to  unity. 

For  laminar  boundary  layer  flow  at  the  blunt-nose  of  a  body  of  revo¬ 
lution  the  heat  transfer  coefficient  can  be  approximated,  for  example,  by 

Slbulkln's  constant  property  formula^’’  applied  to  the  forward  stagnation 

« 

region  behind  the  normal  shock; 


Stj^  -  0.763(G)-’ (pp^Hg)^(Prj^)"0** 


(3) 


If  the  subscripts  X  are  formally  replaced  by  D,  the  corresponding  Stanton 


^Symbols  are  defined  in  Notation 
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number  for  convective  mass  transport  Is  obtained.  This  coefficient  may  then 
be  introduced  into  the  diffusive  contribution  to  the  net  energy  transfer. 

In  order  to  evaluate  the  steady-state  concentration  driving  force  Ac^,  one  must 
invoke  a  knowledge  of  the  "sink  strength"  of  the  surface  for  atoms.  Let  us 
suppose  that  the  atom  recombination  kinetics  at  the  surface  (subscript  w)  are 
described  by  a  first  order  rate  law^  of  the  form»^3,9i ,ea,»4, los 

Theni  in  the  steady  state,  the  conservation  equation  for  atoms  at  the  gas/colid 
interface  (rate  of  consumption  ~  rate  of  supply)  may  be  written: 

(5) 

This  relation  determines  the  "eigen-value"  c.  of  the  atcn:  concentration  esta- 

A  *  W 

bllshed  at  the  wall  if  c,  is  presumed  to  be  known.  If  this  value  for  c. 
is  introduced  into  equation  (2)  and  use  is  made  of  the  fact  that  the  enthalpy 
h  of  the  partially  dissociated  gas  is  comprised  everywhere  of  the  sum^ 
hj.  +  c^Q,  then  the  net  heat  transfer  rate  <1  =  can  be  written  in  the 

form : 

i{  -  0.761  (B0^uJ’^(Pr.^)~°‘^Ah{n-r  (Lej.)°  '  ^-ll^chem}  (6) 

Am 

where  Ah  ,  ~  ♦h  ,  :  $  is  the  correction  factor  C/(HC)j  and  C  is  the 

chem  chem,e 

relevant  catalytic  parameter:’* 

<=  =  “d>‘’  <’'> 


A  reaction  is  said  to  be  first  order  if  this  rate  depends  linearly  on  the 
local  reacta.it  concentration 

With  the  assumption  Q  ~  constant  (see  Section  III) 
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We  have  also  Introduced  the  notation! 

Le,  =  Pr. /Pi'p.  =  "Lewis-Semenov*'^  number"^  (8) 

r  A  u 


It  Is  observed  that  when  C  0,  ♦  0  and  ^  since  the  chemical  enthalpy 

content  h  .  =  c.  Q  of  the  free  stream  is  of  no  consequence  if  the  surface 

chem 1 e  Ate 

is  absolutely  non-catalytlc  (k  -O).  Conversely,  if  there  were  no  atoms  pre- 

w 

sent  in  the  free  stream  (c.  =0),  again  the  heat  flux  <5  reduces  to  4.,  as  it 

A ,  e  A 

would  if  the  heat  of  recombination  Q  were  identically  zero.  In  this  problem, 
three  new  parameters  have  made  their  appearance.  They  are,  respectively! 

(a)  the  catalytic  parameter  C,  defined  by  eq.(7).  (the  ratio 
of  the  characteristic  interfacial  reaction  rate  to  the 
characteristic  convective  diffusion  rate) 

(b)  the  Lewls-Semenov  number  I.e.  -  Pr. /Pr_  =  D/[X/(pc  )] 

[s..  f  \  D  !> 

(c)  t)ie  fraction  of  the  enthalpy  difference  Ah  ~  h^  -  h^  across 

the  boundary  layer  attributable  to  the  chemical  enthalpy 

content,  h  .  -  c,  Q,  of  the  free  stream, 

chem, 0  A ,e 


In  general,  therefore,  each  of  these  parameters  must  be  specified  in  order  to 
calculate  the  heat  transfer  rate.  A  singular  case  arises  when  the  Lewis- 
Semenov  number  LCj,  is  equal  to  unity  since  the  heat  flux  then  becomes  indepen¬ 
dent  of  the  third  of  these  parameters.  It  might  seem,  off  hand,  that  the  heat 
transfer  rate  also  becomes  independent  of  the  chemical  kinetic  (catalytic) 
parameter  C,  but  this  is  not  true  since  the  total  enthalpy  h^  of  the  gas  mixture 
at  the  interface  w  includes  a  chemical  contribution  c,  Q  which  is  a  function 

/\,W 

of  C !  tha t  i  r,  t 


h 

w 


+  n  h 

I  chem,e 


(9) 


The  significance  of  the  subscript  f  will  become  clear  in  discussing  ihe 
opposite  extreme  of  fast  gas  phase  chemical  kinetics 
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where 


1  -  “a.A,.  = 


(10) 


i  - 
r  : 

i-1 


and  we  have  neglected  the  weak  dependence  of  the  chemically  "frozen"  enthalpy 
hj  on  the  diffusion  correction  r|  (see  Section  I'll).  Nevertheless,  it  is  in¬ 
teresting  that  when  Le^  1  the  stagnation  point  heat  flux  q  becomes  explicitly 
proportional  to  the  total  enthalpy  difference  Ah  across  the  boundary  layer  re¬ 
gardless  of  how  much  of  Ah  is  attributable  to  compositional  changes  across  the 
layer.  This  is  often  misinterpreted  as  implying  that  the  heat  transfer  rate 
becomes  independent  of  chemical  kinetic  parameters.  As  discussed  above,  the 
error  consists  In  overlooking  the  fact  that  the  total  enthalpy  h^  of  the  gas 
at  the  Interface  Is  not  known  a  priori,  even  if  the  surface  temperature  is  pre¬ 
scribed.  In  the  case  treated  here,  the  kinetics  of  the  interfacial  atom  recom¬ 
bination  reaction  (together  with  the  surface  temperature)  determines  the  enthalpy 
of  the  gas  at  the  Interface,  and  hence  the  value  of  Ah. 


A  related  misconception  easily  dispelled  is  the  prevalent  notion  that 
the  total  enthalpy  difference  across  the  boundary  layer.  Ah,  is  the  correct 
"driving  force"  for  energy  transport.  Eg. (6)  shows  that  for  the  stagnation 
point,  this  is  only  true  In  the  special  case  l.e^  =  1.  More  generally,  eg. (6) 
reveals  that  the  appropriate  driving  force  is  the  difference  between  a  general¬ 
ized  recovery  enthalpy  and  the  thermal  enthalpy  corresponding  to  surface  tem¬ 
perature.^^’’*  This  can  be  demonstrated  as  follows.  We  first  solve  for  the 
thermal  (frozen)  enthalpy  at  the  surface  which  would  be  reguired  to  cause  the 
total  enei'gy  transfer  rate  4  (in  the  presence  of  chemical  surface  reaction) 

to  vanish.  Setting  ^  =  0  In  eg. (6)  and  solving  for  h.  one  finds* 

f  ,w 


chem 


(11) 


When  the  free  stream  atom  mass  fraction  is  not  small  compared  to  unity, 

qnd  the  thermodynamic  and  transport  properties  of  atoms  and  mole¬ 

cules  a-o  significantly  different  from  one  another,  then  the  heal  transfer 
coefficient  itself  may  couple  appreciably  with  g.  One  aspect  of  this  coupling 
has  been  discussed  recently  by  Inqer  in  reference  (57). 


lU 
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wh*re  .  =  *(0)  h  .  .  This  may  be  considered  to  be  a  recovery  enthalpy, 

chem  chem,e 

with  the  term  ^(0)  (Le^)°**  being  the  effective  recovery  factor  for  free  stream 
chemical  energy.  If  we  now  rewrite  eq.(6)  in  terms  of  this  recovery  enthalpy 

^f,r  ~  ^*^f,w^(!i  =  0  '"®  ^®f  ^ 


4  ^  0.763(pPgPlg)^(Prj^) 


-0.6 


-  "f.w^ 


(12) 


One  rev'ognizes  an  analogy  here  between  the  recovery  of  the  directed 
kinetic  energy  of  the  free  stream  (in  compressible  non-reactlve  heat  transfer) 
and  the  recovery  of  free  stream  chemical  energy  (in  the  present  case).  The 
free  stream  kinetic  energy  at  the  stagnation  point  is  identically  zero,  account¬ 
ing  for  the  absence  of  the  Prandtl  number  Prj^  in  the  (driving  force  for  energy 
transport.  But,  in  general,  both  Pr^  and  Le^  will  appear  in  the  true  driving 
force  for  energy  transfer.  A  revealing  example  is  provided  by  the  flat  plate. 

For  the  case  of  diffusion  controlled  surface  reaction  (<>  1)  the  heat  trans¬ 

fer  distribution  is  found  to  be  given  byj®^*’* 

4  0,332(Re^^^(Prj^)‘f  G  Ah°{l  +  T (Prj^)^-l]  +  [(Le^,)!-!]  — 

Ah  Ah 

2  2 

where  Ah,  ,  “  AC'iii  )  = -i-u  *  and  since  •  1 ,  Ah  .  =  (^c,  =  c,  Q. 

kin  e  chem  A  A,e 

Again,  in  the  special  case  Pr.  ~  1  and  Le,  =  1,  one  could  state  that  the 

difference  Ah°  in  stagnation  enthalpy  is  the  true  driving  force  for  energy 

transport.  However,  more  generally,  the  true  driving  force  is  h^  r  "  w 

wheret 


=  h 


f  .e 


(Prj^) 


h,  ,  + 

kin.e 


(Lef)i 


chem,e 


(14) 


Tlius,  even  when  Pr^^  ^  1  and  Le^  /  1,  the  energy  transfer  distribution  4(x) 
is  given  bys 


4  =  0.332(Re^) 


-i 


(Pr.  )-T  G  (h,  -  h,  ) 

A  f,r  f,w 


(15) 
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Only  when  Prj^  1  and  Le^  1  does  (13)  reduce  to: 

4  =  0.332(Re^)'^(Pr^)‘i  G  (l6) 

Boundary  Lavers  In  Local  Thermocheiiilca  1  Equilibrium 

We  temporarily  leave  the  "chemically  frozen"  case  and  turn  to  the 
opposite  extreme  in  which  the  qas  mixture  within  the  laminar  boundary  layer 
is  everywhere  in  local  therrnochemlca  1  equilibrium.  It  can  be  shown  that  a 
diatomic  gas  In  dissociation  equilibrium  will  behave  as  If  it  were  a  pure 
(single)  substance  with  an  enhanced  thermal  conductivity  (see,  for  example, 
references  15*52,53, 65. 73, 74, 77).  Physically,  the  enhancement  is  the  result 
of  the  diffusion  of  atoms  from  hot  to  cold  regions  of  the  gas  (due  to  the 
change  In  the  equilibrium  atom  concentration  with  temperature)  and  the  subse¬ 
quent  gas  phase  release  of  the  recombination  energy.  If  the  thermal  conduc¬ 
tivity  of  the  equilibrium  mixture  Is  written  \  and  the  chemically  frozen 

^  eq 

thermal  conductivity  is  written  X^,  there  Is  a  simple  relation  between 
X  Af  and  the  corresponding  change  c  /c  ,  in  heat  capacity  attributable 
to  chemical  reaction’^.  In  what  follows,  this  relation  is  used  to  obtain  an 
estimate  of  the  rate  of  heat  transfer  at  the  forward  stagnation  point  of  a 
blunt-nosed  axi -symnit'  L 1 1  c  hooly  when  local  thermocheml ca 1  equilibrium  is 
achieved  everywhere  within  the  gaseous  boundary  layer.  The  Lewis-Semenov 
number  for  atom  diffusion  will  be  assumed  constant*.  This  particular  problem 
is  chosen  because  a  more  rigorous  solution  (in  the  case  of  partially  dissoci¬ 
ated  air)  has  been  obtained^^  by  machine  computation,  so  that  the  accuracy  of 
the  simple  development  given  here**  can  be  checked.  The  present  method, 
furthermore,  provides  3  useful  insight  into  the  way  In  which  rapid  gas  phase 


Computed  as  if  the  composition  did  not  change  with  temperature.  In  view 
of  this  distinction  the  subscript  f  is  implied  on  X  wherever  it  appears  in 
each  of  the  previous  Sections 

This  assumption  hroak?  down  as  the  gas  approaches  the  condition  of  complete 
dissociation,  as  discssed  in  Section  VII 

*Sco,  for  example,  reference  (96) 
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chemical  reaction  should  Influence  the  form  of  heat  transfer  correlation  for¬ 
mulae.  Again*  It  will  be  found  that  the  enthalpy  difference  i&h  across  the 
boundary  layer  Is  the  proper  driving  force  only  In  the  singular  case  of 
le^  =  1.  For  LewlS'Semenov  numbers  different  from  unity*  It  is  Interesting 
to  find  that  the  true  driving  force  for  energy  transport  cannot  be  very  differ¬ 
ent  from  that  obtained  earlier  for  a  catalytic  surface  [eq.(ll)]  In  the  absence 
of  gas  phase  recombination.  This  provides  evidence  In  support  of  a  very  gen¬ 
eral  approach  to  the  problem  of  reacting  boundary  layers*  developed  In 
Section  VIII. 

As  stated  above*  In  a  partially  dissociated  diatomic  gas,  If  thermal 
diffusion  and  other  secondary  diffusion  processes  are  neglected,  the  energy 
flux  vector  can  be  written  In  the  non-reactive  formj 

%  =  -\q  (17) 

where  the  "equilibrium"  thermal  conductivity  X  ^  la  related  to  the  ordinary 

eg  ’ 

frozen  thermal  conductivity  X^  through* 

To  apply  this  result  to  convective  heat  transfer  problems  we  note  that  ordinary 
low  speed  heat  transfer  data  may  be  correlated  in  the  Stanton  form* 


c 


P.avg 


AT 


(19) 


where*  for  gases*  the  dimensionless  Stanton  number  is  usually  represented  as  a 
power  function  of  the  Prandtl  number  Pr^^.  When  dissociation,  diffusion  and  atom 
recombination  occur,  the  changes  in  the  average  specific  heat  and  Prandtl  number 
alone  will  then  tmbody  the  principal  physical  and  chemical  effects.  This  natur¬ 
ally  suggests  the  application  of  eqs.(l8)  and  (19)  but,  for  the  purpose  of  ob¬ 
taining  an  explicit  correlation  pqnatlon,  we  further  introduce  the  average 
properties* 


('^p.eq^avg 


=  Ah/AT 
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(21) 


^°p,f^avg 


=  Ahj/AT 


It  is  through  the  ratio  of  these  two  average  heat 
energy  parameter  of  the  form  explicitly 
when  Le^  ^  1.  The  ratio  of  equilibrium  to  frozen 


capacities  that  a  chemical 
enters  this  class  of  problems 
Prandtl  number  becomes! 


avg 


(22) 


In  applying  this  approach  to  the  axl-symmetric  stagnation  point  heat 
transfer  problem,  use  is  again  made  of  Slbulkln's  laminar  heat  transfer  co¬ 
efficient  [see  eq.(3)],  where,  for  heat  transfer  purposes,  p  p  will  again  be 
introduced  in  place  of  (pM-)^^,^  •  Combining  eqs.  (3),  (l9)  and  (22)  we  then 
Immediately  predict  that  the  equilibrium  heat  transfer  rate  at  the  stagnation 
point  should  be  given  by  a  correlailofi  equation  of  the  form: 


1  •  o  ^ 

^  0.763 Ah  {l  +  (Le^  -  l)-^^}  (23) 

It  will  be  noted  that  the  exponent  (0.6)  on  the  factor  in  brackets,  has  its 
origin  in  the  exponent  (-0.6)  on  the  Prandtl  number  Pr  in  the  non-reactive 

At  t 

heat  transfer  coefficient  [eq.(3)].  For  non-separated  laminar  boundary  layer 
flows  this  exponent  does  not  take  on  a  very  wide  range  of  values  (e.g. .  the 
asymptotic  extremes  0  and  Prj^  ”*00  yield  the  exponents  -l/2  and  -2/3, 

respectively,  corresponding  to  l/2  and  2/3,  respectively,  on  the  bracketed 
"augmentation  factor").  In  general,  we  should,  therefore,  expect  Lewls-Semenov 


^The  computational  results  of  Fay  and  RiddelJ^^  show  that  a  somewhat  better 
choice  would  be: 

(Pli)tvg  =  t(p^)/(p^p^)]°*’ 
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number  factors  of  this  type  to  have  a  weak  dependence  on  the  ratio  of  thermal 
to  vorticity  boundary  layer  thickness.  We  will  return  to  this  point  In  dls” 
cussing  alternate  methods  for  calculating  equilibrium  heat  transfer  rates,  both 
in  laminar  and  turbulent  boundary  layers.  An  alternate  demonstration  of  the 
approximate  validity  of  eq.(22)  can  be  given  by  making  use  of  the  energy  equa¬ 
tion  of  laminar  boundary  layer  theory.  While  this  second  method  is  apparently 
more  restrictive,  it  shows  why  the  problem  may  be  treated  in  terms  of  a  modi¬ 
fied  Prandtl  number. 

The  energy  equation  for  the  laminar  boundary  layer  flow  of  binary 
mixture  of  perfect  gases  may  be  written  in  terms  of  the  static  enthalpy  as 
follows: 


r  ah 


(24) 


Consider  now  the  case  of  a  flat  plate  (dp/dx  ~  0)  with  negligible  viscous  dissl 
patlon  (last  term  of  eq.(24)  small  compared  to  other  terms).  The  right  hand 
side  of  eq.(24)  can  then  be  rewritten  in  tex'ms  ov  enthalpy  gradients  to  obtain: 


(25) 


Inspection  of  this  equation  reveals  that,  if  an  average  value  for  the  quantity 
in  square  brackets  is  introduced,  eq.(25)  reduces  to  the  energy  equation  for 
a  pure  (single)  substance  with  a  modified  Prandtl  number..  A  reasonable  choice 
for  this  "effective"  Prandtl  number  is  seen  to  be: 

"  '■-'f  -  (26) 


where  Ah  ,  is  the  difference  between  h  ,  evaluated  at  the  outer  and  inner 
chem  chem 

edqe  of  the  boundary  layer  and  Ah  is  again  the  difference  h^-h^_  in  the  static 
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enthalpy  (including  chemical  contributions)  across  the  boundary  layer.  Since 
the  heat  transfer  rate  for  a  pure  substance  would  have  been  given  by: 

4  =  St(Re,Prj^)  G  Ah  (27) 

where  St(Re,Prj^)  ~  (Prj^)“^'"^^and  for  most  gases  b  =  •^,  we  conclude  that  for 
comparable  boundary  conditions,  the  heat  flux  In  the  presence  of  equilibrium 
chemical  reaction  will  be  approximately  given  by: 


<S  ~  St(Re,Prj^^P  G  Ah  {]  +(1.8^,  -  (28) 

For  the  case  of  stagnation  flow  this  approach  Is  seen  to  lead  to  the  same 
result  as  obtained  earlier  [eq.(23)].  Graphical  values^  of  the  Lewls-Cemenov 
number  augmentation  factor: 


{l  +  (Le^  -  (29) 

are  shown  in  Fig.  1  for  b  -  j  and  several  values  of  Le^.  A  more  detailed 
discussion  of  the  magnitude  of  the  individual  parameters  and  Ah^j^^^/Ah 
is  postponed  to  n  later  section. 

Comparison  of  eq.(23)  with  eq,(6)  leads  to  the  Interesting  conclusion 

that  if  atom  recombination  does  take  place  the  resulting  heat  transfer  rate  q 

for  a  prescribed  value  of  Ah  and  about  the  same,  regardless  of  whether 

the  recombination  (and  hence  Ah  ,  )  occurs  solely  as  a  result  of  surface  re- 

chem 

action  or  at  equilibrium  within  the  gas  phase.  This  conclusion  is  trivial  In 


On  the  basis  of  comparisons  with  the  computer  solutions  of  reference  35i  the 
procedure  adopted  here  may  be  expected  to  slightly  overestimate  the  Importance 
of  departures  from  the  assumption  Le  =  1 


lfi 
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"?thQ  singular  case  of 
the  factor! 


to  the  factor 


1,  It  la  not  trivial  for  Le^  1.  By  comparing 

{l  t  [(U,)°-‘-  (30) 

{l  +  [(L.,)  -  (31) 


we  find  close  agreement’*  over  a  realistic  range  of  values  of  and 

Lewls-Semenov  number  Le^.  Thus,  heat  transfer  rates  can  bo  estimated  using 
the  driving  force  h,  -  h-  with  h,  ^  given  by  eq.(ll),  regardless  of  the 
location  of  the  atom  recombination  reaction  or  the  magnitude  of  the  Lewis- 
Semenov  number.  So  long  as  the  enthalpy  differences  and  Zih  are  the  same 

In  both  cases,  it  should  not  matter  whether  this  is  the  result  of  gas  phase  or 
surface  reactions,  or  both,  insofar  as  the  heat  transfer  rate  is  concerned. 


An  Alternate  Method  for  Boundary  Lavers  in  Local  Thermochemical  Equilibrium 

Instead  of  using  enthalpy  differences  as  the  starting  point  for  heat 
transfer  calculations,  an  alternate  method  piesents  itself  for  the  case  of 
local  thermo chemical  equilibrium.  This  method  treats  the  prediction  of  heat 
transfer  with  chemical  reaction  as  a  straightforward  variable  property  prob¬ 
lem,  as  done  in  the  previous  section,  but  eliminates  enthalpy  differences  In 
favor  of  differences  in  a  quantity  which  will  be  called,  after  Hansen*^,  the 
"heat  flux  potential".  A  brief  discussion  of  the  use  of  heat  flux  potentials 
(in  place  of  enthalpy  potentials)  for  the  prediction  of  equilibrium  heat  fluxes 
is  given  below. 

In  the  previous  section,  use  has  already  been  made  of  the  fact  that  a 
gas  mixture  in  dissociation  equilibrium  behaves  like  a  pure  substance  with  an 
enhanced  thermal  conductivity  and  specific  heat  which  depend  strongly  on  the 
temperature  and  pressure  level.  Consequently,  the  prediction  of  convective 
heat  transfer  rates  in  such  a  system  has  been  fruitfully  regarded  as  a  variable 
property  problem.  The  groundwork  for  this  approach  was  laid  by  W.  Nernst”  and 


IS 


I.  Larigmuir^^.  Subsequent  developmants  are  aasoclated  with  the  names 
W,  Schotte^°^»  R.  Brokaw"’^®*  C.  F.  Hansen^^*^^,  D,  M.  Mason^***^®  and 
others®*’’*’*’ 

In  clasaical  heat  conduction  theoryi  when  the  thermal  conductivity  X 
is  temperature  dependenti  it  is  convenient  to  define  a  new  dependent  variable 
(p  by  means  of  the  Klrchoff  transformation**’^®: 

r  T 

tp  H  J  X(T)  dT  (32) 

o 

This  dependent  variable  has  appropriately  been  called  the  heat  flux  potential 
since  the  heat  (energy)  flux  vector  Jq  at  any  point  In  the  medium  is  simply  re¬ 
lated  to  the  spatial  gradient  of  q)i 


Jq  =  -  grad  q> 


(33) 


Steady  state  heat  fluxes  in  one-dimensional  stagnant  media  (e»q. «  thermal  con¬ 
ductivity  cells)  are  therefore  directly  proportional  to  the  difference  in  9 
evaluated  across  the  boundaries.  From  the  definition  (32)  It  Is  clear  that, 
for  such  problems,  the  use  of  (p  is  equivalent  to  the  Introduction  of  the  fam¬ 
iliar  temperature  averaged  thermal  conductivity' ’ « 


X 


avg 


(34) 


In  the  case  of  chemically  reacting  gases  we  will  have  X  dependent  upon  pressure 
as  well  as  temperature^  so  that.  In  general,  ip  =cp(T,p).  Since,  at  any  pressure, 
the  thermal  conductivity  X  may  be  regarded  as  the  sum  of  a  chemically  frozen 


Tlie  pressure  dependence  of  the  frozen  thermal  conductivity  Is  small  compared 
to  the  pressure  dependence  of 


iq 


we  may  write*  by  analogy! 


contribution  and  a  reactive  contribution 

f  =  ft  *  <35) 


where 


Xj(T,p)  dT 


(36) 


(37) 


Proceeding  now  to  the  case  of  convective  heat  transfer,  we  recall  that  low 
speed,  nearly  isothermal  (constant  property)  convective  heat  transfer  rates 
in  the  absence  of  chemical  reaction  can  be  correlated  In  the  Nuaselt  formi 


di  =  Nu(Re,Prj^)  X  (Al/l.)  (38) 

Tt  has  been  observed  that  the  principal  effects  of  equilibrium  chemical  re¬ 
action  are  to  change  the  thermal  conductivity  and  specific  heat  of  the  gas 
mixture.  Since  the  thermal  conductivity  appears  explicitly  in  eq.(38),  there 
is  no  question  but  that  the  large  change  in  X  attributable  to  the  dissociation - 
dlf fuslon-qas  phase  recombination  mechanism  must  be  taken  into  account.  In 
contrast,  however,  the  heat  transfer  coefficient  depends  only  slightly 
(fractional  power  law)  on  the  ratio  of  the  specific  heat  to  the  thermal  conduc¬ 
tivity  (via  the  Prandtl  number).  Because  this  dependence  is  usually  weak  to 
begin  with,  and  because  both  the  heat  capacity  and  thermal  conductivity  are  in¬ 
creased  by  comparable  amounts^,  it  has  been  conjectured  (see,  for  example, 
references  65  and  45)  that  in  the  presence  of  chemical  reaction,  the  convective 


If  the  Lewis -Semenov  number  were  identically  unity  X 
would  be  rigorously  equal  to  one  another  in  a  binary 
mixture  (see  eq.flB)) 


a  nH  r  /r 

■p,eq'  ''p,f 
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(39) 


heat  transfer  rate  should  be  given  very  nearly  by 


4  •=  Nu(Re,Prj^^j)(Zi9/L) 


where  the  heat  transfer  coefficient  NUj^  is,  in  this  approximation,  unchanged 
by  the  thermochemical  processes  occurlng  within  the  boundary  layer.  If  so,  this 
would  constitute  an  extremely  potent  computational  technique  since,  apart  from 
the  restriction  to  cases  of  local  thermochemical  equilibrium,  the  method  can 
make  use  of  a  vast  body  of  existing  convective  heat  transfer  data  and  is, 
prima  facie,  free  of  restrictions  as  to  the  behavior  of  the  Lewls-Semenov 
numbers  in  multi -component  gas  mixtures*"^.  One  would  require  only  tabular  or 
graphical  values  of  the  heat  flux  potential  9(1, p),  which  for  a  given  mixture 
could  be  calculated  once  and  for  all,  using  the  methods  outlined  in  references 
13  and  15. 


In  order  to  gain  a  physical  insight  into  the  accuracy  of  this  method, 
as  stated,  as  well  as  to  obtain  an  explicit  estimate  of  the  dependence  of  the 
error  upon  the  known  parameters  and  boundary  data  of  the  problem,  the  following 
procedure’®®  was  adopted.  In  spirit,  the  calculation  closely  parallels  that  of  the 
previous  section. 

If  it  is  conjectured  that  the  heat  transfer  coefficient  in  eq.(39)  is 
to  first  order  Nu(Re,Prj^  ^),  then  this  amounts  to  neglecting  the  effect  of 
thermochemlca 1 ly  induced  changes  in  the  laminar  Prandtl  number  on  the  boundary 
layer  film  conductance.  But,  for  many  boundary  layer  flows,  the  Priindtl  number 
dependence  of  the  non-dimensional  heat  transfer  coefficient  NUj^  is  adequately 

~  (Prj^)*’,  where  typical  values  of  the 
exponent  b  range  from  l/2  to  l/3  depending  upon  the  magnitude  of  the  Prandtl 
number  Pij^  and  the  nature  of  the  flow  field.  Thus,  the  accuracy  of  eq.(39) 
should  be  dependent  upon  the  power  b  of  the  ratio  (Pr.  )  /(Pr.  ,) 

An  explicit  estimate  of  this  factor  in  terms  of  boundary  data  is  possible  by 
introducing  the  quantities: 


represented  by  a  power  law,  i  .e. ,  NUj^ 


(X  „)  =  Affi/AT 

eq  avy 


(40) 
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(41) 


(X.)  =  AmVAT 

f  avg 

If  use  l4  now  made  of  eq,(l8),  together  with  the  assumption  of  constant  Lawis- 
Semenov  number*  we  flndt 


^’’iisslsva  _  fj 


avg 


^cheml 
d/f  J 


(42) 


Here  1*  the  difference  between  evaluated  at  the  outer  edge  of  the 

boundary  layer  (subscript  e)  and  at  the  wall  (subscript  w)  and  Ap  is  simply  the 
difference  in  total  heat  flux  potential  across  the  boundary  layer.  For 

reacting  binary  mixtures  we  conclude  that  the  factor  neglected  when  non- 
reactlve  heat  transfer  coefficients  are  combined  with  heat  flux  potential 
driving  forces  will  he  approximately  given  byi 


{>-['-  E^]  (*3) 

Again,  It  is  Instructive  to  Investigate  an  alternative  derivation  of  this  re- 
sulti  this  time  making  use  of  the  energy  equation  of  laminar  boundary  layer 
theory,  but  expressed  in  terms  of  the  heat  flux  potential  tp  as  the  dependent 
variable.  For  this  purpose  we  make  use  of  a  relation  between  differential 
changes  of  <p  and  h,  readily  derived  from  eg. (18),  1  .e.  i 

[p^]  =  *i{'  ■  ['  -  (4^.) 
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Transforming  the  convective  terms  on  the  left  hand  side  of  eq.(2£i}  with  the 
use  of  eq.(44)«  we  obtaint 


If  an  average  value  for  the  quantity  in  curly  brackets  is  introduced 
eq.(45)  reduces  to  the  heat  conduction  equation  for  a  pure  substance  with  a 
modified  Prandtl  number.  A  reasonable  choice  of  this  effective  Prandtl  number 
is  seen  to  bet 


The  remainder  of  the  argument  leading  to  eq.(43)  parallels  that  given  In  dis¬ 
cussing  the  use  of  enthalpy  as  a  driving  force  for  energy  transport  and  will 
be  omitted  here. 

For  partially^  dissociated  diatomic  gases,  the  Lewis-Semenov  number 
Le^  is  greater  than  unity**,  reflecting  the  fact  that  atom  diffusion  is  a  more 
efficient  energy  transport  mechanism  than  ordinary  conduction  through  the  mix¬ 
ture.  Ttius,  the  factor  (43)  will  be  less  than  unity  for  all  forseeable  cases 
in  which  energy  is  transferred  from  hot  partially  dissociated  gaces  to  cooled 
solids.  Graphical  values  of  this  factor,  for  the  case  b  l/3,  can  be  read 
off  of  Fig. 2.  Inspection  of  (43)  reveals  that,  in  general,  the  error  will 
become  negligible  in  three  distinct  circ  VAHID  WO  nCGS  I 

(a)  if  Le^  -  1 

">>  ‘f  ^  ' 

(c)  If  b  -*  0 


This  assumption  breaks  down  as  the  gas  approaches  the  condition- of  complete 
dissociationj  as  discussed  in  Section  VII 
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Conditions  (a)-and  (b)  have,  parallels  in  evaluating  the  errors  Implicit 
in  the  use  of  total  enthalpy  as  a  driving  force  for  energy  transport  (see  dis¬ 
cussion  of  the  previous  section).  The  magnitude  of  the  Lewis -Semenov  number 
for  atom  diffusion  will  be  discussed  in  Section  VII.  Tlie  magnitude  of 
ZVpchem/^  will  be  discussed  in  Section  III.  devoted  to  the  numerical  importance 
of  Lewis-Semenov  number  "correction"  factors.  It  is  clear  that  as  Le^  “*  1 
the  dependence  of  the  heat  transfer  on  diminishes.  For  any  value  of 

Le^»lf  ^chem^  ^  ^  would  become  proportional  to  the  non¬ 
reactive  driving  force  to  begin  with.  Conversely,  if  ^ 

4  ~  Nu(Re,Prj^  j)  (Le^)  ^  o*"*  equivalently,  d  ^  Nu(Re,Prjj)  where 

Nu(Re,Prjj)  will  be  recognized^  as  the  mass  transfer  coefficient  (Sherwood  num¬ 
ber).  Condition  (c)  will  not  be  encountered  in  the  presence  of  convection. 

It  can  be  concluded  from  this  simplified  analysis  that  the  combination 
of  non-reactive  heat  transfer  coefficients  with  heat  flux  potential  differences 
will  systematically  overestimate  convective  heat  transfer  rates  in  partially 
dissociated  binary  mixtures  by  a  factor  strongly  dependent  on: 

(l)  the  deparLure  of  the  Lewis-Semenov  number,  Le,,  from 
unity  (If  l.e^  >  l)  ^ 

(ll)  the  fraction  of  the  total  heat  flux  potential  differ¬ 
ence  across  the  boundary  layer  attributable  to  chemical 
reaction 

and  weakly  on« 

(ill)  the  conditions  of  convection  within  the  boundary  layer 


It  is  interesting  that,  in  accord  with  eo.(l3),  Pr^  may  be  interpreted  as 
^’^chem  ~  'pjChem  ^^chem 
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A  Formal  Exiorasslon  for  the  Energy  Transport  Driving  Force  In  the  Prasence 
of  Arbitrary  Gas  Phase  and  Interfacial  Reaction  Rataa 


We  have  discussed  the  driving  force  for  energy  transport  In  the  two 
extreme  cases  of  equilibrium  and  chemically  frozen  boundary  layersi  for  Lewls- 
Semenov  numbers  Le^  different  from  unity.  It  has  been  shown  that*  provided 
the  chemical  enthalpy  change  across  the  layer  (due  to  gas  phase  or  surface 
reaction)  Is  the  samei  the  two  effective  driving  forces  are  approximately  equal 
to  one  another.  Encouraged  by  this  resulti  one  is  tempted  to  go  one  step  fur¬ 
ther  and  conjecture  that  the  forgoing  results  are  applicable  even  when  the  re¬ 
actions  in  the  gas  phase  and  at  the  surface^  occur  at  arbitrary  rates.  To  be 
sure*  the  magnitude  of  these  rates  will  Influence  Ah  .  (l.e. .  Ah  .  will 

not.  he  known  a  priori).  Yet  It  would  be  useful  Indeed  if  one  could  statei 
to  a  good  approximation,  homogeneous  and  heterogeneous  reaction  rates  Influence 
the  rate  of  heat  transfer  only  through  the  enthalpy  terms  Ah**,  Ah^j^^i^^/Ah**  (and 

Ahj^ln/^°  when  this  is  not  negligible).  This  amounts  to  neglecting  the  secon¬ 
dary  coupling  effects  which  exist  between  the  equations  of  motion,  species  con¬ 
servation  and  energy. 


Discussion  of  the  prediction  of  (or  the  ^  In  a  multicom¬ 

ponent  gas)  in  terms  of  chemical  kinetic  parameters  will  be  postponed  to 
Section  VI.  We  will  confine  ourselves  here  to  the  conjecture  introduced  above, 
i .e. ,  to  the  development  of  an  approximate,  formal  expression  for  the  energy 
transport  driving  force  in  a  multicomponent  non-equilibrium  system. 


When  the  Lewis  and  Prandtl  numbers  are  Identically  equal  to  unity,  it 
has  been  shown  by  Bromberg  and  I.lpkis’^  and  others*’^  that  this  conjecture  Is 
approximately  true^.  Whan  the  Lewis  and  Prandtl  numbers  are  different  from 


^We  will  consider  here  only  reactions  which  are  catalyzed  by  the  surface}  not 
those  in  which  the  atoms  of  the  surface  are  themselves  reactants 

*For  chemical  reactions  in  the  gas  phase  or  at  the  surface,  and  in  the  presence 
of  mass  addition 
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m: 


11; 


If" 


a  j 

ff-  : 

m..: 


y  •• 

i'  • 


unity*  the  energy  equation  of  lemlnar  boundary  layer  theory  takes  the  formi 


where 


»°5V+Icj[J‘c  dT+h/”'] 


(48) 


This  equation  holda^  in  the  presence  of  chemical  reaction  rates  proceeding 
at  arbitrary  rates  and  in  the  presence  of  free  stream  pressure  gradients. 

For  the  purpose  of  the  following  discussion  the  bracketed  term  {  }■  on  the  right 
hand  side  of  eq.(47)  will  be  rearranged  into  the  formi 


r  , 

P'x*f  ^P%f 


1) 


1) 


^^chem,i 


(49) 


where  the  notation  (  )y  symbolizes  partial  differentiation  with  respect  to  the 
physical  coordinate  y.  Within  the  boundary  layer  the  term  in  square  brackets 
will  vary  from  point  to  point.  With  respect  to  the  establishment  of  the  total 
enthalpy  field  h°(x,y),  it  is  observed  that  this  variation  will  have  an  effect 
similar  to  a  variable  Prandtl  number  except,  in  this  case,  the  average  effec¬ 
tive  Prandtl  number  should  be  the  product  of  some  function  \|/iof  the  para- 

meters  Le^^^,  Ah^hem,!^^”  ^P^,f^avg 


^The  conditions  under  which  the  multi-component  diffusion  terms  can  be  written 
in  this  simple  form  have  been  discussed  by  L.  Lees  in  reference  67 
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Alternativelyi  this  might  be  wrltteni 

=  "’'X.t''’  'SO) 

where,  from  the  form  of  eq.(47).  It  will  be  noted  that  \j/  1  if  Pr.  ,  1  and 

2  * 

Leo  ,  1.  Also,  ilf  1  if  Ah,  ,  0  and  Ah  .  ,  0*  Now  the  function  >1; 

t|l  ^2  kin  cnBiTtfi 

is  not  known  for  the  most  general  case  of  chemically  reacting  compressible 
boundary  layer  flows,  so  that  what  follows  is  to  some  extent  tentative.  Based 
on  the  fact  that  is  known  in  certain  special  cases,  one  can  piece  together 
a  more  general  by  again  making  the  assumption  that  coupling  (interaction) 
effects  are  small  compared  to  the  terms  retainedt 

(1)  For  the  compressible  laminar  flow  of  a  non-reacting  ga8®°**’i 

Ah.  , 

t  •=  1  +  (r,  -  1)  (51) 

*  ^  Ah° 

where  r^^  =  rj^(Prj^  ^  is  the  recovery  factor  for  free  stream 
kinetic  energy. 

(2)  For  the  incompressible  frozen  flow  of  a  fluid  containing 
several  reactants  which  diffuse  toward  a  catalytic  surface 
and  react  there  at  arbitrary  rates,  we  have’**’®*’*i 


Ah. 


,  .  r  /  “"chem,! 

^  V.  •  Ah 


wh 


ere  ’  '’-'f,l)’'‘ 


(52) 

(53) 
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(3)  For  a  binary  gas  in  which  all  the  reaction  occurs  at 
equilibrium  in  the  gas  phassi  it  has  been  shown  that 

the  ijr  given  for  case  (2)  is  approximately  valid. 

z 

(4)  For  the  combined  case  (l)  and  (2)  the  assumption  of  negli¬ 
gible  coupling  leads  to’*« 


^kin  ^  V  f-  ,  \  ^chem.l 


(■■x  -  *  Z  ''D.1  -  *> 


Ah" 


(54) 


fsee,  for  example,  eg. (13)1 

(5)  Since  the  energy  equation  in  terms  of  the  total  enthalpy 
h°  does  not  explicitly  contain  terms  involving  the  kin¬ 
etics  of  gas  phase  reactions,  in  the  most  general  esse  the 
assumption  of  negligible  coupling  should  again  lead  to  a 
value  of  >]/■  not  very  different  from  eq,(54). 

If  this  much  is  accepted  as  plausible,  then  the  energy  eg. (47)  shows  that  the 
bust  choice  of  "driving  force"  for  convective  energy  transport  when 
^  1,  Le^  ^  ^  1  is  \|r  Ah®  or,  approximately: 


i 

Since  Ah®  =  Ah^  +  ^kln  ^  I  ^chem  i’  rewritten: 


Ahj.  + 


r.  Ah,  . 

X  kin 


Ah  . 
chem,l 


(56) 


This  driving  force  ci  n  vanish  (and  hence  the  net  convective  energy  transport 
can  vanish)  under  a  wide  variety  of  conditions.  The  least  interesting  of  these 
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cases  is  that  In  which  each  term  Individually  vanishesi  l.e..  sayi 


kin 


=  0 


^cheiHii”  ® 


(57) 


These  conditions  would  be  satisfiedt  for  example,  when  an  electrically  heated 
non-cata lytic  resistance  thermometer  is  maintained  at  the  flame  temperature 
in  a  nearly  stagnant  mixture  of  combustion  gases.  A  less  trivial  special  case 
of  zero  heat  flux  Is  obtained  by  setting! 

-ihf  =  +  Y,  ro,!  (58) 

i 


For  example,  this  expression  leads  to  the  correct  recovery  temperature  for 
a  catalytic  plate  Immersed  in  a  high  speed  non-equilibrium  dissociated  gas 
stream. 99. '19 

Before  embarking  on  a  discussion  of  the  anticipated  magnitude  of  these 
chemical  effects,  it  should  be  observed  that  in  the  extreme  case  of  local 
thermochemlciil  equilibrium  the  second  law  of  thermodynamics  (or,  qualitatively, 
Le  Chatller'n  principle^)  Imposes  the  condition  that  Ah,  and  Ah  .  should 
have  the  same  algebraic  slgn^.  Thus,  (if  «  o)  an  equality  likei 


~  ^chem 


(59) 


If  a  change  occurs  in  one  of  the  factors  under  which  a  system  is  in  equili¬ 
brium,  the  cystem  will  tend  to  adjust  itself  to  annul,  as  far  as  possible, 
the  effect  of  that  change 

For  nearly  constant  pressure  (transverse)  boundary  layers 
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can  only  hold  If  la  idantlcally  equal  to  zero.  Physicallyi  this  means  tKat 
a  temperature  difference  will  always  ensure  a  net  convective  heat  flux  under 
conditions  of  local  thermochemlcal  equilibrium  (in  low  speed  systems).  This 
Is  consistent  with  the  observation  that  Is  always  a  positive  quantity^. 

In  Section  VIII  we  will  briefly  return  to  this  formal  expression  (56)  for  the 
energy  transport  driving  forcei  and  give  on  alternate  Interpretation  of  Its 
structure*.  We  turn  now  to  an  investigation  of  the  anticipated  effects  of 
non-unity  Lewls-Semenov  number  on  convective  heat  transfer  rates  from  partially 
dissociated  gases  to  solid  surfaces  over  a  range  In  pressures  and  temperatures. 


Tor  dissociation  -1  D^,  (p/TXto)/(RT)3*  c^(l  -  c^),  l.e.i  Is  quad¬ 

ratic  in  the  enthalpy  change  AH  across  the  reaction  Aj  ji.  2A 


Since  the  effect  of  Le^  >  1  la  to  Increase  the  heat  transfer  rate  over  that 
predicted  for  Le^  =  1,  eq.(5)  of  reference  90  should  be  corrected  to  readi 

Nu,  ’C  {l  4  t(u,)P  -  1]  '4^ 

Similarly*  eq.(3-43)  on  pq.45  of  reference  33  should  be  corrected  to  readi 


(Re*)"^  ^  0.70  (Pr*)i  [l  +  (Le"  -  l)~  ] 
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Ill  EFFECTS  DUE  TO  1^6  GREATER  EFFICIENCY  OF  ENERGY 
TRANSPORT  BY  DIFFUSION 

Due  to  its  formal  slmpllcityt  the  singular  case  of  Lewls-Semenov 
number  =  1  has  often  been  used  as  the  basisi  or  starting  point)  for  heat  trans¬ 
fer  calculations.’^^  If  this  is  donei  then  one  must  correct  the  result  for  the 
anticipated  effect  of  departures  from  Le^  =  ll  just  as  one  must  correct  for 
departure  from  ^  ■=  1  when  calculating  heat  transfer  rates  in  the  presence 
of  non-negllglble  viscous  heating.  The  magnitudes  of  these  departures  from 
Le^  =  1  are  discussed  in  what  follows  for  both  equilibrium  and  non-equilibrium 
boundary  layer  heat  transfer. 

It  will  be  recognized  that  the  magnitude  of  these  departures  are  im¬ 
plicit  in  the  results  of  the  previous  sections.  If  we  arbitrarily  decide  to 
use  heat  fluxes  based  on  total  enthalpy  (or  heat  flux  potential)  differences 
as  the  starting  point  Instead  of  using  generalized  "recovery"  enthalpies*®’’*’^*, 
then  a  correction  factor  enters  the  problem  which  is  a  function  primarily  of 
Le^  and  ^®f  '^chem^^  potential  formulation 

la  used  for  non-disslpative  (low  free  stream  Mach  number)  flows.]  In  other 
words,  by  not  using  appropriate  "recovery  enthalpies",  we  forfeit  the  indepen¬ 
dence  of  the  resulting  heat  transfer  "coefficient"  on  such  parameters  as  the 
Lewls-Semenov  number  Le^  and  This  is  analogous  to  the  dependence 

that  the  heat  transfer  "coefficient"  would  have  on  the  Prandtl  number  Pr^^  and 
Ahj^in/Ah°  for  non-reactive  boundary  layer  flows  in  the  presence  of  viscous 
dissipation.  We  will  consider  here  the  expressions  for  the  low  speed  (Mach 
number)  correction  factors  and  discuss; 

(a)  under  what  conditions  the  Lewls-Semenov  number  effects 
should  be  most  noticeable 

(b)  the  anticipated  magnitude  of  these  effects,  based  on 
recent  experimental  data  for  dissociated  air  as  well  as 
combustion  products. 
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Thftt Importanca -  bf  the  lewis -Semenov  Number  Itself 

If  the  enthalpy  difference  Ah  Is  used  as  the  "driving  force"  for  low 
speed  convective  heat  transfer  in  a  binary  reacting  gas  mixture  then  the  appro¬ 
priate  correction  factor  is  seen,  from  eq8.(30)  and  (31),  to  be  approximately! 

{l  t  (r^  -  1)  »  (Uf)’"'’  (60) 

whether  the  boundary  layer  is  In  local  thermochemlcal  eQui librium  or  not. 

Since,  at  most,  Ah  .  /Ah  =  1,  the  maximum  value  this  factor  can  take  on  Is 
chem 

approximately  The  minimum  attainable  value  of  is  zero,  for 

which  the  augmentation  factor  (6o)  becomes  unity  regardless  of  the  magnitude  of 
Le^.  Thus,  in  general! 


1  <  {l  +  [(L.,)'-''  -  1]  }  <  (1.,)’-'“  (61) 


Since  0  <  b  <  1,  it  is  clear  that  the  effect  is  largest  wheni 

(1)  The  Pick  diffusion  coefficient  for  the  energy  containing 
lighter  constituent  is  appreciably  larger  than  the  mean 
thermal  dlffusivity  of  the  gas  mixture  through  which  it 
wanders. 

(2)  A  significant  portion  of  the  enthalpy  difference  across 
the  boundary  layer  is  directly  attributable  to  chemical 
shifts  in  the  composition  of  the  gas  and  not  due  to  tem¬ 
per  a  Uire  (thermal  enthalpy)  differences. 


For  purposes  of  disciisjiion.  (l)  will  be  described  as  a  "favorable  Lewis  number 
condition"*'^.  This  condition  tends  to  exist  for  lean  hydrogen-oxygen  flames, 
for  instance,  since  the  light  hydrogen  atoms  present  wander  through  a  compar¬ 
atively  hejivy  gas  for  which  the  therma 1  dlffusivity  is  correspondingly  reduced. 
In  the  case  of  a  dissociating  diatomic  gas,  the  Lewis-Semenov  number  is  largest 
when  the  atom  concentration  is  smallest,  since  increasing  the  relative  atom 
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conceri,tratlon  hi«  the  tf fact  of  increasing  the  thermal  conductivity  and  hence 
thermai  dlffyalvity  of  the  gaa  mixture.  A  quantitative  example  of  this  behav¬ 
ior  for  the  case  of  pure  hydrogen  will  be  given  in  Section  VII.  When  the  atoms 
are  present  in  "trace"  amounts  at  high  temperaturesi  the  Lewis-Semenov  numbers 
are  always  in  excess  of  unltyt  in  the  case  of  hydrogen*  for  example*  available 
estimates  would  place  it  in  the  neighborhood  of  1.3.  For  weakly  dissociated 
air**  and  nitrogen  tetroxlde” the  Lewis  number  is  estimated  to  be  about 
1.4.  These  values  will  be  nearly  constant  over  a  wide  range  of  temperatures 
in  a  non-equilibrium  (chemically  frozen)  system.  In  an  equilibrium  system* 
however*  this  is  no  longer  true*  since  temperature  and  degree  of  dissociation 
are  no  longer  Independent  (see  Section  Vll). 


Dependence  on  Chemical  Contribution  to  PrivinQ  Force 

It  remains  for  us  to  discuss  the  conditions  under  which  one  would  expect 
^chem^^  to  be  large  (l.e. *  nearest  unity).  For  the  equilibrium  case*  this 
can  be  determined  in  terms  of  temperature  levels  alone,  if  the  total  pressure 
is  prescribed.  The  general  expression  for  If'  dissociating  diatomic 

gas  (a  =  atom*  M  =  molecule)  Isi 


^chem 

Ah 


’’am 

''am*  ^'a 


where  we  have  introduced  the  notatloni 


(62) 


(o)  _  ^(o) 
AM  ■ 


(63) 


Vt,a 


+  (1  - 


T,M 


(64) 


We  first  note  that  when  surface-to-free  stream  temperature  ratio  T  /T  approaches 

w  @ 

zero  and  the  external  temperature  level  is  prescrihed,  Ihen  Ah^^^^/Ah  will 
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approach  the  value  of  In  the  free  stream.  On  the  other  hand,  when 

1  we  have  an  indeterminate  form  for  eq.(62)»  i.e. .  O/O,  which  can  be 
evaluated  using  L'Hospital's  rule: 


Lira  Ah  . 

T  /t  -►  1 

V^e  ^  Ah 


'-is’ 


*  [-lam  *  Xam 


Since  [h^  may  be  identified  with  (Cp  'this  expression 

may  be  written: 


Llm  Ah  .  /Ah  =  [c  .  /c  3  [l  +  (h,  )] 

T  /?  -»  1  chem  '■  p.chem'  p,eq  e  T,AM  AM 

*  tal*  A 


In  the  light  of  these  relations,  consider  a  hypothetical  extreme  case  for  which, 
at  each  pressure  P,  there  exists  a  "threshold"  temperature  T  at  which  the  rel¬ 
ative  atom  concentration  jumps  from  zero  to  unity.  TTien  the  curves  of 
Ah  ,  /Ah,  corresponding  to  the  temperature  ratio  T  /T  “*■  0  and  T  /T  1 

would  show  singular  behavior  when  the  free  stream  temperature  T  passes  through 

ft  ® 

T  (see  Fig, 3b).  In  particular,  when  T  =  T  ~  T  one  would  have 
^  we 

Ah  .  /Ah  -  1|  but  Ah  ,  /Ah  would  be  zero  elsewhere.  On  the  other  hand,  when 
chem  chem 

T  /T  =  0,  Ah  ,  /Ah  would  have  the  behavior  shown  in  Fig. (3b),  jumping  to 

W  0  CnSITl 

some  value  (b^y^pj^/h)  <1  and  then  decreasing  slightly  as  a  result  of  the  almost 
linear  increase  of  h.j.  ^(T)  with  temperature  beyond  T  .  Realistically,  of  course, 
there  will  be  a  threshold  "region"  (interval)  of  temperature  in  which  the  atom 
concentration  will  continuously  vary  from  zero  to  unity  (see  Fig. 3c),  with  the 
resulting  Ah^^g^j^/Ah  behavior  sketched  In  Fig. 3d.  Quantitative  examples  of  this 
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(a)  (c) 


TEMPERATURE,  T  — ►  TEMPERATURE  — ► 

(b)  (d) 


FIG.  3  ATOM  CONCENTRATION  AND  ENTHALPY 
CHANGES  IN  EQUILIBRIUM  DISSOCIATING 
GASES  AT  CONSTANT  PRESSURE 
(a,  b  hypothetical!  c,  d,  actual) 
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^tYlp9  fpr  dissociating  hydrogen  are  given  in  Figs. 4  snd  5.  These  figures  were 
constructed  from  the  data  and  calculations  given  by  Relsfeld^  (reference  83). 
The  behavior  of  ’**^**'  pressure  when  T^T^  =  1  (i.e..  in  the  limit  of 

small  boundary  layer  temperature  differences)  is  shown  in  Fig. 5. 

In  the  case  of  combustion  flamest  the  burned  gas  temperatures  cannot  be 

t 

chosen  at  will  sincei  in  an  adiabatic  system,  the  enthalpy  level  is  determined 
by  the  chemical  energy  content  of  the  metastable  reactant  mixture.  Equilibrium 
energy  changes  in  flame  gases  at  atmospheric  pressure  have  been  discussed  and 
computed  by  Dlxon-Lewls  in  reference  27,  in  particular,  for  the  hydrogen/air 
system  and  for  the  carbon-monoxide/oxygen  system.  Approximate  values  of 
^chem^'^  obtained  from  this  work  show  that  in  the  hydrogen/alr  case  values 
of  for  strongly  cooled  solids  (T^  300°K  -  600°K)  are  In  the  range 

of  10  percent  or  less.  However,  for  solids  maintained  at  temperatures  near 
the  flame  temperature  (2385*^K)  values  of  are  in  the  neighborhood  of 

40  percent.  The  carbon  monoxide/oxygen  flame  gases  exhibit  this  same  increas¬ 
ing  trend  with  Increasing  solid  surface  temperature.  Thus,  heat  transfer  pre¬ 
dictions  based  on  the  assumption  Le^  =  1  (and  using  total  enthalpy  difference 
as  the  driving  force)  should  begin  to  noticeably  underestimate  the  actual  heat 


Owing  to  an  error  in  the  formula  given  In  references  83  and  127  for  calculating 
the  heat  capacity  of  the  reacting  mixture,  we  have  recalculated  the  heat  capa¬ 
cities  tabulated  by  Reisfeld  (reference  83).  The  correct  expression  for  the 
chemical  contribution  to  the  molar  heat  capacity  isi 


'p,chem 


4(1  -ix^)* 


R 


> 


Tabular  values  of  the  heat  capacity,  enthalpy,  thermal  conductivity,  heat  flux 
potential  and  related  functions  for  dissociating  hydrogen  will  be  given  in 
reference  102 


For  the  case 
thalpy  level 
which  can  be 


of  a  high  speed  vehicle  entering  the  earth's  atmosphere,  the  en- 
Is  determined  almost  entirely  by  the  flight  speed  of  the  vehicle, 
arbitrarily  large  (i.e. ,  limited  only  by  the  speed  of  light) 
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i  as  the  solid  8  temperature  approaches  the  flame  gas  temperatures 

ii  ( i  ♦  e « t  a  8  “*■  1). 

j;  Having  discussed  the  conditions  under  which  likely  to  be 

I  large*  it  is  clear  that  similar  considerations  will  apply  to  the  analogous  heat 

■  flux  potential  quantity  For  the  case  of  equilibrium  partially  disso¬ 

ciated  air,  this  has  been  discussed  by  the  writer  in  reference  100.  Here,  use 
was  made  of  shock  tube  recent  data  reported  by  Hansen  in  reference  45.  Hansen 
and  co-workers*^’^*’*’'®  have  displayed  their  experimental  and  theoretical  results 
as  curves  of  q)/cpf  versus  temperature.  In  terms  of  this  quantity,  it  is  readily 
verified  that  for  any  gas  mlxturei 


^^chem 


(t/t,).  -  (\/Tf 


1] 


(68) 


where  it  has  been  assumed  that  the  chemically  frozen  (inert)  contribution 

varies  approximately  as  the  3/2  power  of  the  absolute  temperature.  It  is  again 

of  interest  to  investigate  the  temperature  ratio  extremes  T  /T  0  and 

w  © 

t/t  -*1.  In  the  first  case  by  Inspection  we  obtalnt 

W  0 


w'  ^  e 


(69) 


In  the  second  case  (T  /T  “*■  l),  application  of  L'Hospital's  rule  to  eq.(68) 

W  0 

gives! 


Lim 

T  /T 
w  e 


^chem''^  =  ^  ■  {V'Pf)  [d(cp/(Pf)/d(ln  T)]}’’  (70) 
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Of*  Is  the  fact  that  these  two  limits  become  ehual  only  where 

d(^/<j)j)/dT  =  0^.  This  Is,  of  coureo,  true  for  temperature  levels  at  which 
there  Is  no  dissociation!  but  it  can  also  be  true  at  higher  values  of  the 
temperature,  l.e.i  In  regions  of  nearly  complete  dissociation.  For  Inter¬ 
mediate  temperatures  T  dissociation  causes  d(q>/cp^)/dT  >  0  and  hence,  In  this 
range t 


Llm 

V’', 


w'  e 


(71) 


An  inranedlate  consequence  of  the  foregoing  reasoning  when  applied  to  the  case 
of  air  Is  the  property  that.  If  the  free  stream  stagnation  temperature  level 
la  In  the  neighborhood  of  nearly  complete  oxygen  dissociation,  l.e..  If 
dC(p/<p^)/dT  =  0,  then  the  Lewla-Semenov  number  correction  factor  will  be  the 
same  In  either  extreme  T  /t  0  or  T  /t  “*'1  and  given  approximately  byi 


where  Le^  will  have  a  value  near  1.4-.  For  intermediate  values  of  the  boundary 

layer  temperature  ratio  T  /t  ,  this  Is  no  longer  true. 

w  © 


t 


since: 


d(<p/<Pf) 

d(ln  T) 


and  X^l/p^  >  0  for  T  >  0,  it  follows  that  this  derivative  vanishes  when 

®/cp.  =X  ^Ar  “  ^  Le^r  (c.  /c  ,)  -  l]  (for  a  dissociating  diatomic  gas) 
eq  T  r  ■  p,eq  p,r 
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Using  ‘R®l8fi3,ci'!i  calculatiani**  for  pure  hydrogen  as  a  starting  p^lnt^, 
we  have  constructed  ®*  *  function  of  temperature  and  pressure  level 

for  the  two  extreme  cases  =0»1.  Hiese  are  shown,  respectively,  in 

^  Figs.  6  and  7.  Tabular  values  of  the  input  data  will  be  given  in  a  separate 

report’®*.  Computations  for  the  case  of  multi -component  combustion  gases  would 
f  lead  to  qualitatively  similar  results  but  would  considerably  be  more  tedious 

and  time  consuming  to  carry  out.  Nevertheless,  with  the  advent  of  electronic 
computation  machines,  and  the  theoretical  methods  outlined  in  references  15  and 
13,  this  is  now  well  within  thw  realm  of  possibility.  As  usual,  one  must  weigh 
+he  magnitude  of  the  effort  l...olved  and  its  ultimate  utility  against  the  un¬ 
certainties  in  the  values  of  much  of  the  Input  date.  For  the  more  common  gases, 
there  is  no  question  in  the  writer's  mind  but  that  these  calculations,  if  pro¬ 
perly  displayed,  would  be  of  great  value  in  making  accurate  heat  transfer  pre¬ 
dictions  in  high  temperature  systems  of  future  interest. 

;  Paralleling  the  reasoning  given  earlier,  it  can  be  concluded  that  con- 

I 

vective  heat  transfer  predictions  based  on  the  assumption  Le^  =  1  (using  heat 
flux  potential  driving  forces  with  non-reactlve  heat  transfer  coefficients) 
should  tend  to  noticeably  overestimate  the  actual  heat  fluxes  as  the  solid  sur¬ 
face  temperature  T  approaches  the  flame  gas  temperature  T  (i .e. ,  as  T  /T  l). 

W  6  yy|  ^ 

The  anticipated  magnitude  of  this  effect  can  be  read  of  Fig.  2  If  one  estimates 
Uj  .nd 

For  the  case  of  non-eouillbrlum  boundary  layers,  longer 

a  function  of  temperature  levels  alone,  since  the  kinetics  of  both  the  gas  phase 

and  surface  recombination  reactions  will  determine  the  value  of  Ah  .  /Ah  es- 

chem 

tablished  at  each  point  along  a  solid  surface.  As  an  example,  consider  an 
electrically  heated  platinum  resistance  thermometer  being  maintained  at  the 
flame  gas  temperature  in  a  rich,  hydrogen-oxygen  flame  at  low  pressures.  Avail¬ 
able  experimental  evidence  shows  that,  within  and  immediately  behind  the  primary 
reaction  zone,  the  hydrogen  atom  concentration  is  larger  than  that  corresponding 


See  footnote,  page  30 
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to  local  thermo chemical  equilibrium  by  one  or  more  orders  of  magnitude,*®  In 
a  situation  such  as  this  would  be  unity  so  long  as  there  is  any  yas 

phase  or  surface  hydrogen  atom  recombination!  This  example  also  suggests  that 
while  the  Lawis-Semenov  number  conditions  themselves  may  be  more  favorable  in 
lean  flames*’’',  the  Lewls-Semenov  number  augmentation  factors  can  be  appreciable 
in  rich  flame  gaaes  as  well  because  of  the  larger  values  of  Ah^^gi^/Ah,  Cl.e., 
despita  the  fact  that  Le^  is  reduced).  An  interesting  feature  of  the  limiting 
case  Ah^j^gj^/Ah  ~  1  is  the  fact  that  the  heat  flux  becomes  proportional  to 
St(Re,Pr|^)  where  St(Re,Prjj)  will  be  recognized  as  the  mass  transfer 

coefficient  (Stanton  number  for  convective  diffusion).  Physically,  this  corres¬ 
ponds  to  energy  transfer  by  diffusion-surface  reaction  mechanism  alonej  i  .e. , 
in  the  absence  of  ordinary  convection.  The  fraction  Ah^j^^^/Ah  as  well  as  the 
enthalpy  difference  Ah  Itself  will  he  determined  in  part  by  chemical  kinetic 
and  aerodynamic  factors.  In  the  case  where  no  gas  phase  recombination  occurs, 
then  Ah^^g^/Ah  will  depend  on  a  catalytic  parameter  similar  to  that  discussed 
earlier.  As  an  example,  we  recall  the  hypersonic  stagnation  point  problem  and 
note  thatt 


^chem 

Ah 


♦h 


ch  em ,  e 


“’t  *  •''ch.m,. 


where 


♦  =  C/(l  +  C) 


h 


chem,e 


c 


A,e 


(73) 


(74) 

(75) 


Here,  we  have  again  assumed  that  the  surface  reaction  obeys  first  order  kin¬ 
etics,  When  gas  phase  recombination  cannot  be  neglected  ^  itself 


The  weak  dependence  of  Ah^  on  C,  which  would  be  caused  by  compositional  changes 
if  the  heat  capacities  of  the  atoms  ai>d  mulecules  were  suf ficien''ly  different, 
is  neglected 
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will  depend  in  a  more  complex  manner  on  both  gas  phase  and  surface  chemical 
kinetic  parameters.  Since  each  of  these  parameters  contains  "aerodynamic” 
variables  as  well,  we  have  here  a  situation  in  which  the  dependence  of  the 
heat  transfer  rate  on  a  change  in,  say,  mass  velocity  G  =  p  u  is  no  longer 
confined  to  the  heat  transfer  coefficient  itself,  but  extends  into  the  esta¬ 
blishment  of  the  heat  transfer  driving  forces  themselves.  A  discussion  of 
this  dependence,  however,  will  be  postponed  to  Section  VI. 


45 


IV  CALCULATION  OF  THE  TURBULEWT  FILM  CONDUCTANCE 


IN  AXI-SYMMETRIC  NOZZLES 


As  suggested  by  the  discussion  of  Section  II,  It  Is  useful  to  divide 
the  convective  heat  transfer  problem  into  two  parts  (the  determination  of  a 
film  conductance  and  a  driving  force)  since  it  is  readily  demonstrated  that 
the  film  conductance  (or  heat  transfer  coefficient)  has  a  relatively  weak  de¬ 
pendence  on  the  details  of  the  physical  property  value  profilesn  hence,  rela¬ 
tively  crude  techniques  for  taking  these  property  variations  into  account  are 
often  quite  adequate  from  a  practical  point  of  view.  As  a  first  approximation, 
it  la  therefore  reasonable  to  assume  that  the  heat  transfer  coefficients  will 
be  substantially  the  same  as  those  determined  experimentally  or  semi-theoreti- 
cally  for  the  case  of  non-reactive  boundary  layer  flows.  In  extreme  cases  the 
accuracy  of  this  approach  may  be  questionable}  yet,  it  is  to  be  expected  that 
the  dominant  functional  relationships  cannot  be  very  different  from  those  ob¬ 
tained  using  this  technique.  Thus,  one  has  a  simple  starting  point  to  which 
refinements  can  be  added  when  necessary. 

In  the  previous  sections  wa  have  discussed  that  part  of  the  problem 
which  1*.  (in  the  case  of  rapid  gas  phase  or  surface  reactions)  xelatlvely  in¬ 
sensitive  to  fluid  mechanical  (aerodynamic)  conditions.  We  turn  now  to  a  brief 
consideration  of  the  calculation  of  film  conductances  with  emphasis  on  turbulent 
boundary  layer  development  in  axl-symmetric  rocket  motor  nozzles.  In  contrast 
to  the  determination  of  energy  transport  driving  forces  when  the  chemical  kin¬ 
etics  are  rapid,  determination  of  the  film  conductance  is  primarily  an  aero¬ 
dynamic  problem. 

Fig.  8  shows,  schematically,  the  meridional  cross-section  of  a  conven¬ 
tional,  regenera tively  ccolei  liquid  propellant  thrust  chamber.  The  reactants 
"A"  and  "B"  are  supplied  at  constant  rates  to  the  chamber}  however,  reactant  A 
Is  first  made  to  pass  through  cooling  passages  which  envelop  the  nozzle  and 
combustion  chamber.  Each  reactant  enters  the  combustion  chamber  itself  through 
an  array  of  injectors  (atomizers)  at  station  i.  By  the  time  the  nozzle  inlet 
(station  c)  is  reached,  vaporization,  mixing,  and  chemical  energy  release  are 
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assumed  to  be  nearly  complete.  The  hot  reaction  products  are  then  expanded  to 
ambient  pressure  in  a  converging  -  diverging  nozzle  (effuser)  having  a  physical 
throat  section  (station  t)  at  which  the  mass  velocity  G  =  pu  necessarily  passes 
through  a  maximum.  To  prevent  "burnouts"  due  to  melting,  oxidation  or  erosion, 
the  combustion  chambv;r  and  nozzle  walls  must  be  maintained  at  temperatures 
(usually  less  than  1200'\)  which  are  appreciably  smaller  than  the  combustion 
qas  temperatures  (usually  2500\  -  4000°K).  This  implies  that  the  cooling 
system  must  be  able  to  cope  everywhere  with  the  resulting  gas-side  heat  flux. 

For  long  thrust  durations  in  motors  which  are  sufficiently  large,  efficient 
cooling  systems  can  be  designed  which  make  use  of  the  heat  capacity  of  one  of 
the  propellants  ,  (usually  the  fuel)  before  it  enters  the  chamber.  Ideally, 
this  can  be  accomplished  without  having  to  admit  any  of  A  into  the  regions  of 
the  expansion  section  where  the  combustion  gas  pressure  is  lower  than  the  cham¬ 
ber  pressure  P^.  This  then  corresponds  to  the  familiar  case  of  heat  transfer 
to  solid  surfaces  which  are  internally  cooled  but  exposed  externally  to  high 
temperature,  high  velocity  gases.  Tiie  remainder  of  the  present  discussion  will 
be  directed  solely  at  the  gas-side  heat  transfer  problem  under  these  cooling 
conditions.  Accurate  methods  for  predicting  the  gas-side  heat,  transfer  co¬ 
efficient  could  then  be  used  as  the  starting  point  for  an  overall  nozzle  cool¬ 
ing  design  study,  such  as  that  performed  by  Curren,  Price  and  Douglass**,  for 
the  case  of  high-performance  chemical  rockets,  and  by  Robbins,  Bachkin  and 
Medeiros**  for  the  case  of  nuclear  rockets. 

While  it  is  generally  agreed  that  the  prediction  of  heat  transfer  rates 
in  the  combustion  chamber  itself  has  been  impeded  by  gross  uncertainties  in  the 
flow  pattern  caused  by  the  propellant  injection  and  heat  release*',  in  the  words 
of  Bartz’^'i  "somewhere  between  the  entrance  to  the  nozzle  and  the  nozzle  throat 
it  Is  expected  that  convection  due  to  average  one-dimenciona 1  flow  will  begin 
to  dominate  the  problem.  Under  these  conditions  Lheie  is  hope  for  success  of 
analytical  predictions  of  heat  flux  based  only  on  considerations  of  convection 
due  to  the  average  one-diniensiona !  gns  flnw".  Since  this  is  precisely  wnere 
the  cooling  problem  becomes  most  critical  (near  the  throat)  and  since  the  ex¬ 
pansion  (diverging)  section  constitutes  a  large  fraction  of  the  total  nozzle 
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surface  requiring  cooling,  the  problem  is  not  an  academic  one. 

The  experimental  and  theoretical  determination  of  local  heat  transfer 
coefficients  in  rocket  nozzles  has  occupied  numerous  investigators  for  the  past 
20  years.*’**’ With  the  emergence  of  less  conventional  nozzle  designs  for 
space  flight  applications,  there  has  been  renewed  effort  in  the  search  for  a 
simple  but  yet  sufficiently  accurate  correlation  formula‘s*  to  describe  the 
distribution  of  heat  transfer  coefficient  in  terms  of  both  geometric  properties 
of  the  nozzles  and  fluid  properties.  For  this  purpose,  Spalding''^  and  Mayer"^* 
have  Independently  made  use  of  Ambrok's  approximate  boundary  layer  analysis*  to 
arrive  at  an  expression  for  the  local  Stanton  number  which  combines  the  attri¬ 
butes  of  reasonable  accuracy  with  simplicity.  In  general,  it  is  to  be  recom¬ 
mended  over  the  use  of  modified  pipe  flow  formulas^ (based  on  local  nozzle 
diameter)  and  boundary  layer  methods  based  on  the  Reynold's  analogy*’*®  (in 
which  the  skin  friction  distribution  is  determined  prior  to  the  application  of 
Reynold's  analogy  between  skin  friction  and  heat  transfer).  The  former  method 
Is  known  to  underestimate  the  heat  flux  in  accelerating  flows  (favorable  pressure 
gradient)  while  the  latter  method  overestimates  the  heat  flux  in  accelerating 
flows.  By  combining  the  Ambrok  method*  with  the  Rubesin-Eckert  reference  tem¬ 
perature  method**’*^’**  (to  correct  for  the  effects  of  variable  fluid  pronerties) 
Mayer’’*  has  developed  the  following  relation  for  the  local  Stanton  number 
Stj^(x)  in  an  axl -symmetric  nozzle  of  local  diameter  d(x)j 

-I  3  -3  X  *  1 

Stj^(x)  =  0.0296(Prj^^  j)  d^O*^  [G(x)/p;]dx}  ’ 

0 

t 

where 

0  ^  ’  (p^/pg)  ’  ~  3  +  '‘O 


If  the  ahsn’i'it.e  visnosity  |i  of  the  nas  mixt.\ire  varies  as  and  the  density  p 
varies  as  1/T  .  - 

*The  viscosity  coefficient  appearing  in  the  integrand  in  eg. (76)  is  to  be 
evaluated  at  the  outer  edge  of  the  boundary  layer 


(76) 


(77) 


Kr 


(78) 


G(x)  = 

and  the  subscript  *  corresponds  to  evaluation  at  a  reference  temperature 
which  in  the  present  case  should  be  given  implicitly  byi 

h^(T^j)  =  0.50[hf(T^)  +  hj(T^)]  +  0‘22[rj^(Prj^^  (79) 

where  rj^(Pr^  ^)  is  the  recovery  factor  for  free  stream  kinetic  energy.  For 
laminar  boundary  layers  on  blunt  nosed  hypersonic  vehicles,  Eckert  and  Tewflk*'^ 
have  recently  shown  that  the  reference  enthalpy  concept  coupled  with  Lees 
analysis**  for  constant  specific  heat  c^  and  density-viscosity  product  pji 
predicts  local  convective  heat  transfer  coefficients  which  are  in  satisfactory 
agreement  with  the  results  of  more  exact  boundary  layer  calculations.  In  the 
present  case,  It  is  therefore  tentatively  suggested  that  the  heat  transfer 
coefficient  St^(x)  given  by  eg. (76)  be  combined  with  the  local  enthalpy  differ¬ 
ence  Ah°  and  corrected  for  effects  of  non-unity  Lewis  and  Prandtl  number  by  the 
factor  I 

where 


^kin 


(81) 


chem 


(82) 


Here  we  have  made  use  of  the  following  approximate  values  (see,  for  example, 
reference  28)  for  the  recovery  factors  and  r^  for  free  stream  kinetic  energy 
and  chemical  energy,  respectively J 

1 

^^^X  f^^  (turbulent  flow)  (83) 


SO 


(laminar  or  tu-rbuient  flow) 


(84) 


In  cases  where  the  diffusion  of  more  than  one  light  "chemical  energy  carrier" 
is  to* be  taken  Into  account*  it  is  tentatively  proposed  that  the  i?«t  term  in 
eq.(80)  be  replaced  by  the  sumt 

2  Ah 

y  [(Le.  y  -  1]  (85) 

^  Ah° 

in  accord  with  the  discussion  of  Section  II.  This  expression  should  be  formally 
valid  for  arbitrary  values  of  both  the  gas  phase  and  surface  chemical  kinetic 
parameters,  alnce  these  parameters  primarily  establish  only  the  magnitude  of 
Ah°  and  the  This  will  be  discussed  further  in  Section  VI.  We  first 

turn  to  the  simplest  case  of  local  thermochemical  equilibrium  in  the  gas  phase. 
For  this  purpose  one  must  have  available  extensive  thermodynamic  charts  for 
the  propellant  system  in  question. 
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FIG.  9  ENTHALPY  VERSUS  MIXTURE  RATIO  CHART 
FOR  THE  PRODUCTS  OF  THE  REACTION  OF 
A  AND  EL 
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V  THERMODYWMIC  CALCULATION  OF  ENTHALPY/MIXTURE-RATIO  CHARTS 


Engineering  calculations  for  systems  in  which  "real"  fluid  effects  are 
appreciable  are  facilitated  by  the  availability  of  thermodynamic  data  in 
graphical  form.  Unfortunately,  there  is  no  universal  representation  of  thermo¬ 
dynamic  data  which  is  equally  convenient  for  all  calculations.  In  the  case  of 
heat  transfer  c jlculations,  two  types  of  charts  have  found  widespread  use: 

1.  enthalpy  versus  entropy  at  constant  mixture-ratio 

2.  enthalpy  versus  mixture  ratio  at  constant  pressure. 

Tne  first  is  the  well-known  Mo' Her  diagram,  most  useful  for  systems  of  fixed 
overall  clienilcal  (atomic)  composition.  On  the  other  hand,  in  systems  where 
two  individual  reactants  may  be  combined  in  arbitrary  proportions^,  it  has  been 
pointed  out  by  Lutz*’,  Reichert**  and  more  recently  hy  flpalding  and  co-work- 
gf3li3,l14,l13  that  entha Ipy/mixture-ratlo  diagrams  (at  fixed  total  pressure) 
are  often  more  useful.  A  particularly  interesting  discussion  of  the  appli¬ 
cation  of  this  type  of  diagram  to  problems  of  heat  and  mass  traiifer  in  aero¬ 
nautical  engineering  Is  given  in  reference  115. 

The  basic  features  of  such  a  presentation  pertinent  to  the  present 
discussion  are  displayed  schematically  in  Fig. 9.  Tlie  extremetles  on  the 
abscissa  correspond  to  pure  "A"  and  pure  "B",  W^en  only  two  types  of  atoms 
(or  nuclell)  are  present  In  the  system,  only  one  chart  is  necessary  and  this 
scale  can  be  chosen  In  such  a  way  that  the  composition  parameter  measures  the 
fractional  composition  (by  mass)  of  one  of  the  atoms  to  the  total,  regardless 
of  the  particular  chemical  state  of  aggregation.  This  choice  gives  rise  to  a 


Design  studies  for  regeneratively  cooled  engines**  suggest  that  under  some 
marginal  cooling  conditions  mixture-ratio  compromises  may  be  necessary.  For 
a  "fuel-cooled"  engine,  this  would  generally  call  for  a  richer  mixture  than 
that  giving  maximum  specific  impulse 
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chart  with  the  following  properties t 


1.  The  state  point  e  corresponding  to  steady  flov;  combustion  of 
A  and  B  at  the  total  pressure  in  question  can  be  located  from 
the  Intersection  of  the  straight  line  joining  the  input  en¬ 
thalpies  h^»  hg  and  the  vertical  line  giving  the  mixture  ratio 


2.  All  states  within  a  gaseous  boundary  layer  will  lie  along  the 
line  of  constant  mixture  ratio,  since  extra-nuclear  chemical 
reaction  within  the  layer  will  only  change  the  chemical  state 
of  aggregation  and  not  the  overall  atomic  composition  of  the 
mixture.  In  particular,  the  state  point  W  corresponding  to 
any  prescribed  wall  temperature  can  be  identified  by  loca¬ 
ting  the  isotherm  for  T^.  The  stagnation  enthalpy  difference 
Ah°  between  the  combustion  product  in  the  free  stream  and  at 
the  wall  is  then  read  off  the  ordinate  directly. 

3.  The  adiabatic  flame  temperature,  when  needed,  can  be  obtained 
by  identifying  the  Isotherm  passing  through  point  e.  It  is 
significant,  however,  that  this  need  not  be  done  to  determine 
the  driving  force  for  energy  transport. 

If  the  Isotherm  corresponding  to  the  surface  temperature  T  is  not  in  the  region 

w 

of  dissociation  for  the  lowest  pressures  encountered  in  a  rocket  motor,  then 
to  a  good  approximation,  the  enthalpy  difference  determined  in  the  above 
manner  would  apply  everywhere^  within  the  thrust  chamber  (at  all  pressure 
levels)  since  the  "free  stream"  stagnation  enthalpy  h°  does  not  change  during 
the  iso-energetic  expancion  process  In  the  noz/'.le  (despite  the  fact  that  the 
chemical  rearrangement  and  temperature  changes  occur).  However,  when  the  chem¬ 
ically  frozen  Prandtl  and  Lewis  numbers  are  different  from  unity,  additional 


^For  a  nearly  Isothermal  wall 
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information  is  needed  to  calculate  local  heat  transfer  rates,  since  the  way  in 

which  the  energy  is  partitioned  (directed  kinetic  energy  as  well  as  chemical 

and  thermal  energy)  across  the  boundary  layer  is  no  longer  immaterial.  This 

is  displayed  by  the  appearance  of  the  terms  (r^^  -  l)  Ahj^j^^/Ah°  and 

(r_  .  -  l)  Ah  .  in  the  semi-theoretic  correlation  obtained  by  combining 

D  y  1  cn  6in  f  1 

eq.(76)  with  (30) .  To  properly  take  these  terms  into  account,  one  would  have 
to  supplement  entha Ipy/mlxture-ratio  charts  of  the  type  described  above  with 
Molller  diagrams  and  equilibrium  composition  diagrams.  Unfortunately,  there 
is  neither  enough  experimental  data  available  to  justify  taking  account  of  these 
terms  nor  to  Justify  their  neglect.  The  discussion  of  the  previous  sections 
may  therefore  be  useful  in  shedding  some  light  on  the  question  of  when  these 
additional  terms  should  become  important.  This  writer  ventures  the  quess  that 
the  effort  may  he  justified  when  dealing  with  nozzle  liners  (ceramic  inserts, 
coatings,  etc.)  used  in  conjunction  with  high  performance  propellant  combina¬ 
tions,  since  these  techniques  tend  to  increase  the  interior  surface  temperatures 
in  the  nozzle  and  thereby  increase  the  fractions  ^kin^'^’° 

appearing  in  eg.(80). 


Entha Ipv/Mlxtur 9-Ratio  Charts  for  Hvdrogen/Oxv 
of  lU,  30  and  6u  Atmospheres 


len  Combustion  at  Total  Pressures'. 


Using  the  thermochemical  data  compiled  by  Huff,  Gordon  and  Morrell®* 
(1951),  we  have  constructed  large  enthalpy/mixture-ratio  diagrams  for  oxy- 
hydroqen  combustion  at  three  total  pressures!  10,  30  and  60  atmospheres  (see 
foldouts  1,  2  and  3).  The  mixture-ratio  parameter  chosen  for  the  ab¬ 

scissa  is  defined  as^! 


weight  (or  mass)  of  oxygen  (any  form)  in  mixture 
total  weight  (or  mass)  of  mixture 


In  Appendix  I,  the  relation  between  this  parameter  and  alternate  mixture-ratio 
parameters,  in  current  use,  is  given. 
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and  enthalpies  per  unit  mass  are  given  on  the  ordinate  in  the  units  kilocalories 
per  gram.  To  avoid  negative  values  of  the  absolute  enthalpies  (i .e. .  thermal  + 
chemical)^,  the  convention  of  reference  56  has  been  adopted  by  arbitrarily 
assigning  the  absolute  enthalpy  of  zero  to  H2O  and  Oj  in  the  solid  (crysta l) 
phase  at  zero  degrees  Kelvin.  This  will  have  no  affect  on  the  graphical  calcu¬ 
lation  of  enthalpy  differences.  The  remaining  requisite  thermochemica 1  data  is 
given  in  Table  1  shown  below  and  in  reference  56. 


TABLE  1 

Assignment  of  Absolute  Enthalpies  (Kcal/Mole)  In  the  Oxy-Hvdroqen  System 


Substance 

Formula 

Phase 

TempC^K) 

Ah, 

form 

h 

Hydrogen 

Hz 

Gas 

0 

67.4169 

Hydrogen 

Ha 

Ga  s 

293.16 

0 

69.4407 

Hydrogen 

H 

Ga  5 

0 

51.62 

85.3285 

Oxygen 

Oa 

Ga  s 

298.16 

0 

4.1109 

Oxygen 

Oa 

Ga  s 

0 

2.0362 

Oxygen 

Oa 

Crysta 1 

0 

2.0362 

0.000 

Oxygon 

0 

Ga  8 

0 

53.586 

59.6041 

Wa  ter 

HjO 

Gas 

293.16 

-57.7979 

13.6938 

Wa  ter 

HjO 

Gas 

0 

11.3311 

Wa  ter 

H,0 

Crysta 1 

0 

-68.4350 

0.000 

Hydroxyl 

OH 

Gas 

0 

10.0 

4^1.7266 

The  "sum  of  the  sensible  (frozen)  and  chemical  contributions  to  the  enthalpy 
has  been  called  the  "reaction  enthalpy"  in  references  49»  69  and  32.  In  the 
light  of  current  practice,  it  is  recommended  that  this  terminology  be  dropped, 
ot.horwlsp  specified,  enthalpy  should  always  be  taken  to  include  the 

chemical  contributionj  1 . e, ,  h  -  h.,  H  h  , 

-  f  chem 
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Computations  were  performed  using  an  Iterative  method  on  an  IBM  650  digital  com¬ 
puter.  The  equilibrium  equations  were  left  in  their  non-linear  form.  Only  six 
major  species  were  considered  to  exist  in  the  mixture  of  combustion  gases j  i  .e. . 

H  Ha  0  Oa  HaO  OH 

and  their  Individual  physical  and  thermal  properties  were  assumed  to  combine  in 
accord  with  the  laws  governing  mixtures  of  perfect  gases.  Departures  for  non- 
ideality  were  not  taken  into  account  for  any  of  the  pressures  or  isotherms  shown. 
Also  Included  on  these  plots  are  two-phase  regions  with  a  limited  number  of 
isotherms  . 

While  more  recent  thermochemical  data  has  since  become  available  (see, 
for  example,  reference  50),  and  further  refinements  can  be  incorporated  in  the 
calculations  in  the  future  (particularly  at  the  higher  pressure  levels),  the 
graphical  results  are  probably  sufficiently  accurate  for  convective  boat  trans¬ 
fer  calculations  in  chemical  propulsion  systems  utilizing  the  combustion  of 
hydrogen  and  oxygen.  A  series  of  sample  calculations  making  use  of  these  charts 
will  be  given  in  a  later  report.  A  summary  of  useful  supplementary  calculations 
on  the  hydrogen-oxygen  system,  which  may  be  found  in  the  current  literature,  is 
given  in  Appendix  1. 

As  in  the  case  of  Mollior  (enthalpy-entropy)  diagrams  for  reacting  mix¬ 
tures,  the  detailed  chemical  composition,  while  needed  to  construct  the  diagram, 
does  not  make  Its  explicit  appearance  on  the  diagram  itself.  Tnus,  it  Is  not 
possible  to  use  the  chart  alone  to  determine,  say,  across  the  boundary 

layer,  due  to  gas  phase  hydrogen  atom  reassociatlon.  If  the  Lewis-Semenov  num¬ 
ber  is  expected  to  be  greatly  different  from  unity  (see  Section  VIl),  it  has 
already  been  pointed  out  that  this  additional  information  will  be  of  interest. 


Isotherms  in  the  two  phase  regions  are  straight  lines,  since  a  point  moving 
along  such  an  isotherm  represents  a  change  only  in  the  relative  amounts  of  the 
two  states  points  at  the  extremities  of  the  isotherm 


57 


In  these  cases,  enthalpy- (atomic)  composition  data  should  be  supplemented  by 
actual  species  composition  data.  This  is  also  needed  for  purposes  of  computing 
transport  properties  (see  Section  VIl)  and  total  gas  density  (or  molecular 
weight).  Changes  in  mixture-ratio,  pressure  and  temperature  level  cause  dram¬ 
atic  changes  in  the  chemical  composition  of  the  gases  in  the  boundary  layer, 

as  can  be  seen  from  Figs. 10,  II  and  12.  Fig. 13  shows  the  mole  fractions  x„  and 

H 

’'HjO  hydrogen  atoms  and  water  vapor,  respectively,  for  a  stoichiometric 
mixture  (f  =  0.838)  over  a  range  in  pressure  level.  Tabulated  chemical  compo¬ 
sitions  and  mean  molecular  weight  for  the  stoichiometric  case  are  included  in 
Appendix  2  at  temperatures  between  1000°K  and  4000°K  (200*^K  intervals)  for  each 
of  the  total  pressures  10,  30  and  60  atmospheres.  A  brief  discussion  of  the 
computation  of  the  transport  properties  of  such  gas  mixtures  will  be  postponed 
to  Section  VII.  In  this  connection,  it  should  be  observed  here  that  at 
T  =  3400*^K,  p  =  10  atm.  (say)  molecular  fragments  account  for  some  18  particles 
per  hundred,  with  hydroxyl  radicals  (OH)  accounting  for  10  of  them^.  This  in¬ 
troduces  an  additional  uncertainty  into  estimates  of  the  high  temperature  trans* 
port  properties’**,  since  the  requisite  cross-sections  are  not  well  known  for 
collisions  involving  free  radicals.  It  is  hoped  that  this  uncertainty  will  be 
reduced  as  a  result  of  a  sustained  experimental  and  theoretical  effort  in  this 
direction. 


See  Table  6,  Appendix  2,  for  actual  species  composition  data  in  the  combustion 
of  a  stoichiometric  mixture  at  10  atmospheres  total  pressure 


mole  fractkm 


FIG.  10  PRODUCT  GAS  COMPOSITION  VERSUS  MIXTURE  RATIO  PARAMETER 
FOR  OXYHYDROGEN  COMBUSTION;  P  =  IO  ATM. 


FIG.  II  PRODUCT  GAS  COMPOSITION  VERSUS  MIXTURE  RATIO  PARAMETER 
FOR  OXYHYDROGEN  COMBUSTION;  P  =  30  ATM. 
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VI  EFFECTS  OF  CHEMICAL  NON-EQUILIBRIUM  WITHIN 


THE  FREE  STREAM  AND  WITHIN  THE  BOUNDARY  LAYER 
ON  CONVECTIVE  HEAT  TRANSFER  IN'  ROCKET  MOTORS 

A  good  deal  of  attention  has  been  directed  in  the  past' to  the  effect  of 
chemical  non-equilibrium  on  the  specific  impulse  (thrust  per  unit  mass  flow)  of 
rocket  engines.  In  the  inviscid  (core)  flow,  this  chemical  non-equilibrium  may 
be  the  result  of; 

1.  high  cooling  rates’^'^’^®  (large  rate  of  change  of  temperature 
with  time)  on  the  vicinity  of  the  nozzle  entrance  section  com¬ 
bined  with  large  rates  of  density  decrease,  causing  dissociated 
combustion  products,  initially  at  equilibrium,-  to  lag  behind  in 
attempting  to  follow  these  aerothermodynamic  changes 

2.  failure  to  achieve  thermochemical  equilibrium  in  the  combustion 
chamber  itself,  prior  to  the  expansion  process®’  (e.q. ,  in¬ 
complete  combustion). 

Relatively  little  attention  has  been  paid  to  the  anticipated  effects  of  non- 
equilibrium  on  convective  heat  transfer  rates.  In  this  case,  compositional  lags 
cm  occur  within  the  boundary  layers  themselves  as  well  as  in  the  core  (free 
stream)  flow.  The  effects  of  chemical  non-equilibrium  and  the  conditions  under 
which  these  effects  should  be  most  marked  will  be  briefly  discussed  in  the  pre¬ 
sent  section.  Here  again,  it  is  useful  to  keep  in  mind  recent  studies  of  non¬ 
equilibrium  flows  over  blunt-nosed  hypersonic  bodies  (in  particular,  see  ref¬ 
erences  18,  19  and  40) .  In  qualitatively  comparing  these  f lows^’’^*'^®,  the  stag¬ 
nation  point  on  a  vehicle  can  be  likened  to  the  nozzle  entrance  section,  the 
sonic  point  to  the  throat  section,  and  the  after-body  to  the  divergent  expansion 
section.  While  the  analogy  is  not  complete  in  every  detail,  the  points  of 
similarity  are  frequent  enough  to  make  the  comparison  illuminating. 


Chemical  non-equl librium  can  alter  the  heat  transfer  rates  everywherci 
primarily  through « 

1.  the  resulting  Increase  of  enthalpy  of  the  gas  at  the  interface, 
and  hence,  the  reduction  of  the  enthalpy  difference  Ah°  across 
the  boundary  layer 

2.  the  resulting  changes  of  the  j/Ab® 

3.  the  resulting  change  in  the  transport  properties  of  the  mixture 
(and,  hence,  the  film  coefficient)  caused  by  alterations  in  the 
temperature  and  chemical  composition  fields 

4.  the  resulting  change  of  Ah.  ,  /Ah®. 

Kin  • 

Of  these,  it  is  expected  that  the  first  will  be  the  most  important.  Wo  have 
seen  that  while  the  second  and  fourth  "mechanisms"  disappear  for  those  cases 
in  which  the  Lewis  and  Prandtl  numbers  are  near  unity,  the  first  mechanism  re¬ 
mains  due  tn  the  "damming  up"  of  energetic  species  along  the  interface.  This 
will  occur,  however,  only  when  neither  the  kinetics  of  either  the  gas  phase  or 
Interfaclal  reassoclation  reactions  are  able  to  cope  with  the  capacity  of  the 
free  stream  to  supply  these  species. 

From  a  practical  point  of  view,  the  most  important  implication  of  the 
lag  in  gas  phase  reassoclation  rates  is  the  possibility  that  the  heat  transfer 
rate  will  become  sensitive  to  the  nature  of  the  heat  transfei-  surfaces  (see, 
for  example,  reference  25) »  not  because  of  surface  temperature  or  emlssivlty, 
but  because  each  material  will  have  a  different  ability  to  catalyze  the  recom¬ 
bination  of  molecular  fragments  incident  upon  it. 

Chemically  Fru/tin  Boundary  Lavers  with  Catalytic  Surface  Reaction 

To  an  accuracy  which  is  probably  consistent  with  the  application  of 
eg. (76)  to  the  rocket  motor  problem,  it  is  possible  to  make  a  quantitative  es¬ 
timate  of  the  effect  of  finite  surface  activity  for  the  extreme  case  in  which 
no  reassociation  occurs,  either  in  the  free  stream  or  within  the  boundary  layer. 
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The  method  used  here  corresponds  to  the  Heymann  and  Frank-Kamenetzkll  quaslsta- 
tionary  method^  of  classical  heterogeneous-diffusional  kinetics.****'^’’-'' 

Let  Yj  the  recombination  probability’**  for  the  labile  species  1  in¬ 

cident  upon  the  wall  w  In  question.  We  then  define  a  set  of  "reaction  vslo- 

citles"  k  ,  by^’j 
w,i  ' 


-  t(kTj/(27nn^)]  Yj 


(87) 


and  assume  that  the  equilibrium  concentration  of  species  1  at  the  wall  temper¬ 
ature  Is  neqliqlble.  The  mass  balance  equation  for  species  i  may  be  written  asj 


St„  ,  G  Ac,  =  k  .pc, 
D,i  1  w,l‘^w  l,w 


(83) 


where 


Stjj  ,  local  mass  transport  coefficient  (Stanton  number) 

*  for  species 

G  “  local  mass  velocity  pu  In  the  free  stream 


In  terms  of  a  set  of  new  streamwise  coordinates’**^®! 


-1 


(89) 


we  then  have  the  following  relation  for  the  quantities  Ac,/c.  i 

1  1 « e 


(90) 


Considered  in  the  light  of  recent  applications  of  boundary  layer  theory  to 
bodies  of  arbitrary  shape,  this  Is  equivalent  to  a  "local  similarity"  assump' 
tlon  (see,  for  examole,  references  19,  66,  60) 
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From  low 
obtained 

with  the  appropriate  ratio  of  Reynold's  analogy  factors}  In  this  case  by 

(Le.  ,)*.  Thus,  the  heat  transfer  distribution  would  still  be  given  by  the 
*  1  * 

following  equation! 


speed  heat-mass  transport  similitude  theory,  each  of  the  Stj^  ^  can  be 
from  the  heat  transfer  coefficient  St^(x)  [eq.(76)3  by  multiplication 


Ah, 


a(x)  =  Stj^(x)  G  Ah°  {l  +  (rj^  -  l)-~  + 


Ah'- 


D,i 


j^^^chem.i 
Ah‘ 


•}  (91) 


In  this  case,  howeven 


^^chern,!  ^1  ^'^chem,i^e 


(92) 


Ah  ^Ahj.  +Ah^^^  + 


I'l'" 


chcm,l^e 


(93) 


It  will  be  noted  that  If  each  of  the  Interfaclal  reaction  velocities  k  , 

i 

approached  infinity,  then  each  of  the  coordinates  Zj  would  become  infinite  re¬ 
gardless  of  the  physical  distance  x  downstream.  In  this  extreme,  by  eq.(90), 
each  of  the  catalytic  activity  corrections  would  approach  unity  and,  as  a 
result,  the  enthalpy  changes  Aii®  and  the  Ah  .  ,  would  take  on  their  maximum 

CnGm 1 1 

values.  The  energy  transfer  rate  in  this  extreme  should  not  be  very  different 

from  that  corresponding  to  the  opposite  extreme  of  equilibrium  recombination 

within  the  boundary  layer.  However,  for  any  realistic  material,  eq.(87)  shows 

that  the  reaction  velocities  k  ,  cannot  be  infinite  and,  depending  on  the 

iiiaqniLude  of  the  recombination  probabilities  y, ,  each  k  ,  can  take  on  all 

values  between  0  (non-cata lytic)  and  [kT  /(2rTm, )]*  (perfectly  catalytic).  When 

each  of  the  reaction  velocities  k„  ,  is  zero.  eg. (91)  would  still  applyi  however, 

1 

since  =  0,  we  havet 


Ah  .  ,  =  0 

v#n  6in  f  1 


(94) 
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(95) 


Ah°  =  Ahf  +  Ah 


kin 


This  gives  the  minimum  local  heat  flux  everywhere  andi  hence,  the  minimum  total 
(integrated)  heat  flux.  Between  these  two  extremes,  the  heat  flux  can  take  on 
all  intermediate  values  depending  on  the  values  of  the  reaction  velocities  k 

w,i 

for  each  nozzle  material  chosen.  It  should  be  noted  that  the  streamwise  co¬ 
ordinates  Zj^(x)  depend  upon  aerodynamic  factors  as  well  as  chemical  kinetic 
factors.  Thus,  an  increase  in  chamber  pressure  level  has  the  same  qualitative 
effect  as  an  increase  in  the  chemical  activity  of  the  surface  (roughly  speaking, 
the  coordinates  z^  will  depend  on  the  chamber  pressure  raised  to  the  l/5  power). 
For  similar  reasons,  the  effect  of  increased  physical  scale  is  also  to  increase 
the  coordinates  Zj^.  Each  of  the  latter  changes,  however,  decreases  the  llkll- 
hood  that  no  reassoclatlon  will  occur  in  the  gas  phase. 


F.ffact  of  Gas  Phase  Chemical  Kinetics  on  Heat  Transfer  to  Non-Cata lytic  Surfaces 

The  extreme  sensitivity  of  the  convective  heat  flux  to  the  catalytic 
activity  of  the  nozzle  surfaces  indicated  by  the  foregoing  analysis  will  not  be 
realized  in  practice  if  homogeneous  (qas  phase)  recombination  reactions  in  the 
bounds  I y  layers  themselves  are  not  negligible.  In  the  limit  of  very  fast  homo¬ 
geneous  kinetics,  this  sensitivity,  in  fact,  virtually  disappears^.  However, 
before  qualitatively  discussing  the  interaction  of  both  gas  phase  and  surface 
recombination  rate  parameters,  we  turn  to  the  effect  of  finite  qas  phase  recom¬ 
bination  parameter^  in  the  absence  of  surface  catalysis.  Even  here,  however. 


For  arbitrary  surface  temperature  and  chemical  surface  activity,  the  equili¬ 
brium  diffusion  fluxes  will  not  be  compatible  with  the  kinetics  at  the  inter¬ 
face.  This  gives  rise  to  a  thin  "non-equilibrium  sub-layer"  within  the 
generally  near-equilibrium  boundary  layer,  but  does  not  appreciably  Influence 
the  overall  heat  flux’* 

*For  a  multi -component  qas  mixture,  there  will  be  more  than  one  gas  phase  chem¬ 
ical  reaction  parameter  even  if  certain  equilibria  are  established  among  sev¬ 
eral  of  the  species.  For  the  purposes  of  discussion,  however,  we  will  consider 
a  single  gas  phase  chemical  rate  parameter 
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our  discussion  must  at  prossnt  remain  qualitative. 

Just  as  the  influence  of  arbitrary  interfacial  kinetics  is  indicated  by 
the  magnitude  of  a  non-dimensional  parameter  In  which  chemical  kinetic  and  aero¬ 
dynamic  factors  are  combined,  the  Influence  of  gas  phase  kinetics  will  be  in¬ 
dicated  by  the  magnitude  of  a  non-dimensional  parameter  which,  likewise,  contains 
chemical  kinetic  apd  aerodynamic  factors.  Physically,  each  of  these  parameters 
may  be  interpreted  as  the  ratio  of  a  chemical  rate  to  a  diffusion  rate  or, 
equivalently i,  the  ratio  of  a  diffusion  time  to  a  chemical  (life)time, 

Thus,  pressure  level  and  physical  scale  again  make  their  veiled  appearance. 

Fig. 14  attempts  to  qualitatively  show  how  the  local  heat  transfer  rates  to 
non-catalvtic  thrust  chamber  walls  are  apt  to  vary  with  chamber  pressure  level. 
The  ordinate  shown  is  a  measure  of  the  fractional  "recovery"  of  chemical  energy 
in  the  free  stream  by  virtue  of  reassociation  reactions  within  the  developing 

boundary  layers  ^r.hem^^chem  eq  ^®f  ^  ^*^® 

pressure  Is  sufficiently  large,  departures  fiom  Lh'*  equilibrium  heat  transfer 
rate  should  be  negligible  throughout  the  chamber.  For  smaller  chamber  pressures 
departures  will  begin  to  occur,  noticeably  in  the  region  of  the  nozzle  throat, 
and  tend  to  persist  further  downstream.  At  still  smaller  pressures,  these  de¬ 
partures  should  become  noticeable  eveiywlieie.  Finally,  we  have  the  extreme  in 

which  the  gas  phase  reactlcns  are  "frozen"  (p  “*  0)  causing  Ah  .  /Ah  ,  _  to 

c  cnem  cnem,eq 

vanish  everywhere. 

At  present  it  is,  unfortunately,  not  possible  to  Insert  accurate  numerical 
values  on  a  set  of  curves  of  the  type  sketched  in  Fig. 14-  There  are  cases  cited 
in  the  literature  In  which  non-equilibrium  effects  (specificity  to  surface 
material)  are  apparently  observed  under  conditions  for  which  one  would  have  ex¬ 
pected  large  values  of  the  gas  phase  recombination  rate  parameter^.  On  the  other 
hand,  recent  heat  transfer  Ha+a  for  flat  plates  exposed  to  the  flow  of  oxy- 
acetylene  flame  gases  at  one  atmosphere  total  pressure  suggest  that  the  molecular 


^If  the  free  stream  is  not  in  thermochemical  equi librium  Ah 
*^c.heiTi,e  ^chem,eq,w 


chem,eq 


implies 
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FIG.  14  NORMALIZED  CHEMICAL  CONTRIBUTION  TO  THE  CONVECTIVE  HEAT 

FLUX  TO  A  NONCATALYTIC  WALL 


fragments  present  in  the  external  flow  (h  atom  mole  fractions  of  up  to  approxi¬ 
mately  10  percent)  do  not  survive  their  journey  through  the  turbulent  boundary 
layer5’»'°^.  If  anything  is  to  be  concluded,  it  is  that  additional  investigation 
on  all  fronts  is  needed  in  order  to  build  up  a  body  of  experience  (however 
idealized)  adequate  for  making  quantitative  estimates  of  these  when  they  are 
apt  to  be  important.  Having  discussed  the  extremes,  there  remains  the  inter¬ 
mediate  case  in  which  the  gas  phase  and  surface  recombination  rate  parameters 
toasther  detsp.lne  If  Is  Itself  not  nsgllglbl. 

then  the  energy  transfer  becomes  potentially  sensitive  to  both  parameters. 

Several  qualitative  features  of  this  type  of  problem  can  be  extracted  from  the 
theoretical  investigations  of  Grelflnger^,  and  Hirschfelder^* ’ 109  fgj. 
case  of  a  conductivity  cell,  of  Chung’*  for  a  Couette  flow  model,  and  Scala’°* 
in  the  case  stagnation  point  heat  transfer.  Apart  from  the  conclusions  already 
discussed,  namelyt 


1.  the  heat  transfer  becomes  insensitive  to  the  qas  phase  re¬ 
combination  parameter  when  the  surface  recombination  parameter 
becomes  very  largej 

2.  the  heat  transfer  rate  becomes  insensitive  to  the  surface  re¬ 
combination  parameter  when  the  gas  phase  recombination  parameter 
becomes  very  largei 

3.  the  heat  transfer  rates  In  cases  1  and  2  are  very  nearly  equal» 


one  can  anticipate  some  "coupling"  effect  for  Intermediate  values  of  each  para¬ 
meter,  since  the  catalytic  surface  reaction  infhiences  the  qas  phase  reaction 
by  causing  a  reactant  depletion  "at  a  distance".  It  further  appears  from  the 
work  of  reference  I6  that  when  both  parameters  take  on  intermediate  (comparable) 
values,  the  sensitivity  of  the  heat  flux  to  changes  in  the  surface  recombination 
parameter  is  greater  than  the  sensitivity  to  changes  in  the  gas  phase  recom¬ 
bination  parameter.  Also  of  interest  is  the  anticipated  importance  of  the  sur¬ 
face  temperature  level  both  in  influencing  the  local  rates  of  gas  phase 
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reactions^®  (mainly  by  changing  in  local  density),  as  well  as  in  Influencing 
the  magnitude  of  the  surface  recombination  probability®^.  The  present  dis¬ 
cussion  suggests  to  the  writer,  at  least,  that  if  workable  approximate  corre¬ 
lation  formulae  arc  within  the  realm  of  possibility,  they  will  be  of  the  form 
of  eq.(9l)  except  that  the  functions  introduced  in  our  treatment  of  the 
chemically  frozen  case  will  become,  more  generally: 


®l(zi 


(G) 


y  M  (w) 

z  1  , z*  , 


(w)  ) 

i  ’  •  •  •  ^ 
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that  i:5,  functions  of  both  gas  phase  (G)  as  well  as  surface  (w)  chemical  re¬ 
action  rate  parameters.  For  a  binary  dissociated  qas  mixture,  we  already  know 
something  about  from  the  work  referenced  above.  Thus: 


4  (  OO  ,  Z^'^^)  ~  1 

(93) 

(99) 

«(0,  zM)  t  z/^^) 

(100) 

<1(7.  ,  0)  as  in  I  nferences  3f)  and  40 

(101) 

If  the  coupling  effects  described  above  were  small,  one  might  expect  sviper- 
position  formulae  of  the  following  type  to  be  approximately  valid: 

•  ~  1  -  [l  -  ♦(z^^\o)]  [1  -  ♦(0,z^”^)]  (97a) 


In  view  of  the  limited  size  of  the  existing  dictionary  of  solutions,  interim 
mcasureb  of  Lliif'  type  nitr  not  unjustified.  It.  is  i.oped,  however,  that  more  ef¬ 
fort  and  ingenuity  will,  in  tlie  future,  go  into  tlie  correlation  of  computer 
solutions  so  that  tae  spot  calculations  which  one  so  f:-'pgnent ly  encounters  in 
the  literature  (i .e. ,  for  "typical  re-entry  conditions")  can  be  tied  together 
into  an  intelligible  pattern. 
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VII  ESTIMATION  OF  THE  LEWIS-SEMENOV  NUMBER  AND  OTHER  SIGNIFICAOT 
MOLECULAR  TRANSPORT  PROPERTIES 


Since  chemical  reaction  rates  are  explicit  functions  of  reactant  con¬ 
centrations  and  the  temperature,  these  rates  will  be  implicit  functions  of  all 
molecular  parameters  which  Influence  the  establishment  of  the  concentration  and 
temperature  fields  within  the  boundary  layer.  This  dependence  on  molecular 
transport  properties  exists  even  in  the  case  of  turbulent  boundary  layers,  since 
it  is  known  that  the  major  resistance  to  the  transport  of  heat  and  mass  occurs 
In  a  relatively  thin  sub-layer  near  the  wall  in  which  molecular  transport 
effects  are  dominant.  Across  this  sub-layer,  temperature  and  species  concen¬ 
trations  undergo  their  largest  variations.  In  this  section,  we  turn  our  atten¬ 
tion  to  the  estimation  of  the  molecular  transport  properties  in  high  temperature 
gas  mixtures,  with  emphasis  on  the  determination  of  the  Lewis-Semenov  number. 

For  a  more  detailed  discussion  of  the  estimation  of  viscosity  and  Prandtl  num¬ 
ber  Pr.  ,,  as  well  as  the  corresponding  accuracy  requirements,  the  reader  is 

A,  t 

referred  to  earlier  papers  in  this  field  (references  7  and  8). 

Viscosity  and  Chemically  "F rozen"  Heat  Conductivity 

For  pure  gases,  viscosities  can  ordinarily  he  measured  with  greater  ease 
and  accuracy  than  heat,  conductivities.  In  most  cases,  these  data  can  be  extra¬ 
polated  over  a  wide  range  of  elevated  temperatures  using  Chapman-Enskog  theory*’ 
combined  with  a  realistic  Interaction  potential.  Assuming  then  that  each  of 
the  component  viscosities  p^(T)  is  known,  in  an  N-component  mixture  of  composi¬ 
tion  Xi,X2,  ...  x^,  ...  Xj^j  (>nole  fractions),  a  number  of  formulae  are  avail¬ 
able  for  estimating  the  effective  viscosity  of  the  mixture.  The  reader  has  no 
doubt  already  observed  that  these  vary  widely  in  their  convenience  and  rigor. 

A  stmi-empirlca 1  formula  which  combines  reasonable  accuracy  with  computational 


69 


convenience  is  that  of  Wilke^**t 


>* '  I  [  i 


i=i 


j=i 


(102) 


where 


=  (2)'^  [1  + 


(nijAj)] 


1 

■T 


1  1 
[i  +  (ii^/'nj)^(nij/ni, )  ] 


(103) 


As  stated  above,  thermal  conductivities  of  the  Individual  species  at  high 
temperatures  are  probably  most  , accurately  obtained  from  the  corresponding  vis¬ 
cosities  In  conjunction  with  a  kinetic  theory  law  of  the  Eucken^*  type^  i 

^1  ^  (Cv,iVli/^i)t(i)(7.032Y  -  1.720)]  (104) 


Here  v  is  the  ratio  C,  ,/E  ,  of  molar  heat  capacities*.  Again,  for  combining 

1  V  1 1 

the  constituent  thermal  conductivities  ,  one  has  available  many  "mixing  rules 
of  varying  accuracy  and  convenience,  Tn  the  writer's  opinion,  a  rational  for¬ 
mula  which  combines  reasonable  accuracy  with  computational  convenience  is  that 
recently  derived  by  Mason  and  Saxena^^j 


II 


X 


r 


F.  {',^[F/j^i]'’} 


1^1  j-i 


(105) 


This  type  of  relation  is  most  successful  for  diatomic  molecules.  For  poly¬ 
atomic  molecules  (e.o. ,  water  vapor)  experimental  conductivity  data  should  be 
used  when  possible 

*AccurrtLe  determination  of  llie  individual  heat  capacities  are  possible  by  the 
application  of  statistical  mecnanics  to  available  spec Lroscopi c . da ta  (see,  for 
example,  reference  tio) 
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where  the  are  Identical’  to  those  qlven  above  [eq.(l03)].  The  formal  sim- 
llarlty  of  eqs.(l02)  and  (lOO)  makes  this  method  particularly  useful. 

Plf fusion  Coefficients  for  Molecular  Fragments 

Due  to  the  complexity  of  the  rigorous  laws  of  multi -component  (N  >  2) 
diffusion^"*’ it  is  felt  that  the  present  state  of  knowledge  does  not  warrant 
their  attempted  use  for  problem*  of  the  type  described  herein.  Fortunately,  in 
many  casoc,  the  molecular  fragments  which  are  the  "thermochemical  energy 
carriers”  can  be  regarded  as  present  In  "trace"  amounts.  When  this  Is  a  reason¬ 
able  approximation,  the  effective  Pick  coefficient  for  the  diffusion  of  the 
labile  species  i  through  the  mixture  Is  given  byt 

N  -1 

where  the  are  the  bi nary  diffusion  coefficients.  In  a  two-component  system 
^1  mix  binary  diffusion  coefficient  D)j,  It  can  be  shown,  more¬ 

over,  that  in  this  special  case  diffusion  coefflcl(»nt  is  virtually  independent 
of  concentration. 

While  many  moderately  high  temperature  binary  diffusion  coefficients  are 
accurately  known  for  stable  gas  pa  1 r8’*°* ’**»  this  is  not  yet  true  when  one  of 
the  components  is  a  molecular  fragment.  Even  in  the  case  of  the  most  widely 
studied  of  these  coefficients  (i .e. .  that  pertaining  to  the  diffusion  of  hydro¬ 
gen  atoms  through  diatomic  hydrogen),  it  will  be  seen  that  considerable  uncer¬ 
tainty  remains. 

In  1928,  Harteck'f*  dei;ermined  the  vlaccsity  of  H/Hj  mixtures  of  known 
composition^  by  applying  Poiselle's  equation  to  the  flow  of  dissociated  mixtures 

^Through  the  use  of  a  molecular  effusion  (Wrede-Harteck)  gage  technique 
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through  a  small  tube.  Amdur^  (1936)  has  examliieU  these  data  and  calculated  the 
corresponding  binary  diffusion  coefficient  on  the  assumption  that  both 

H  and  Ha  behave  as  van  dei  Waal’s  gases  (hard  elastic  spheres  of  constant  dia¬ 
meter  with  attractive  forces  varying  as  an  inverse  power  of  the  distance)^. 

'Allien  hydrogen  atoms  are  present  in  trace  amounts,  Amdur  gives  the  relationt 


=  3.335  X  10 


(tVp)  [1  +  (31.9/t)] 


1 


(107) 


where  T  is  the  absolute  temperature  in  K,  p  is  the  total  pressure  in  atmos¬ 
pheres  and  the  resulting  coefficient  has  the  units  cm* (sec)"’.  Although 

Harteck's  data  apply  only  from  273°K  to  373°K,  Amdur's  estimate  is  considered 
to  be  reasonably  reliable  to  temperatures  of  about  600°K. 


Prior  to,  as  well  as  following,  the  Harteck-Amdur  estimate  of  D„  ..  , 

n  -rlj 

there  have  been  several  different  values  calculated  for  this  coefficient.  Three 
such  estimates  are  compared  to  eq.{l07)  in  Fig. 15  over  the  temperature  range 
300°K  to  4000°K.  The  curves  marked  a,b,c,  were  each  calculated  using  the 
Lennard-Jones  6tl2  nonpolar-nonpolar  interaction  potential,  but  with  different 
size  and  energy  parameleis.  These  parameters  are  Mctod  in  Table  2. 


TABLE  2 

Assumed  Lennard-Jones  6il2  Interaction  Potentia 1  Parameters 
for  Atom-Molecule  Diffusion  in  Hydrogen  Cl  =  a  torn. 2  =  molecules) 


Author 

Reference 

ffia (a) 

Eia/K  ('’k) 

L.  Lees  (1956) 

66 

2.637 

38.0 

R.  Brokaw  (l96o) 

12 

2.798 

38.0 

M.J.  Roisfeld  (1957) 

83 

2.94 

30.79 

?.  B'jther'anH  pntpn+lal 
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It  is  observed  that  Lees'  choice  of  interaction  potential  parameters^  would  re¬ 
sult  in  the  largest  values  of  D„  „  as  well  as  the  steepest  temperature  deoen- 

H-Ha 

dencBj  whereas  the  extrapolation  of  the  Harteck-Amdur  relation  leads  to  the 
lowest  diffusion  coefficients  and  the  weakest  temperature  dependence.  In  1912, 
Langmuir*^  made  the  a  priori  estimate j 


1  -1 

„  =0.514  X  lO’*  (tVp)  [1  +  (77/T)]  (108) 

At  temperatures  above  lOOo'^K,  this  leads  to  diffusion  coefficients  which  fall 
reasonably  close  to  those  given  by  Brokaw  (b)  and  Relsfeld  (c).  Recent  measure¬ 
ments  by  Wise  and  co-workers'*^*’*^'  lead  to  high  temperature  diffusion  coeffi¬ 
cients*  which  fall  very  close  to  the  extrapolated  Harteck-Amdur  curve  (d). 

An  alternate  way  of  displaying  the  spread  In  the  available  estimates  of 

D„  „  is  to  compare  the  ogulvalent  hard-sphere  diameters  of  the  hydrogen  atom 
H'Hj  0 

at  300  K.  Several  results  collected  at  random  from  the  literature  are  given  in 
Table  3. 


Based  on  the  estimate  that  CTta  should  be  nearly  equal  to  the  difference  between 

CTaa  and  throe-eighths  of  the  equilibrium  Internuclear  separation  distance  of 

the  parent  dlatotnlr.  molecule;  and  ^ajA 

t  ^ 

Corresponding  to  a  hard  snhere  effective  diameter  Ou  of  2.43  A 

(US.,  1) 
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TABLE  3 


Assumed  "Hard  Sphere*'  Hydrogen  Atom  Diameters  (T  =  300^K) 


ty(A) 

Reference 

1.9 

Lees**,  L.  (1956) 

2.2 

Brokaw^*,  R.  (i960) 

2.3 

Reisfeld**,  M.  (1^57) 

2.5 

Harteck^*,  P.  (192B)  Amdur^,  I.  (1936) 

2.1 

Bonhoeffer,  K.F.,  Harteck,  P.  (1933)^° 

1.1 

Tanford’’*,  C.  (1947) 

2.4 

Semenov’’®,  N.N.  (1958) 

2.2 

Wise’***’* 5,  H.  (1959) 

2.0 

lanqmulr*^,  I.  (1912) 

1.9 

Warren’*’,  D.R.  (1952) 

These  were  computed  by  findinq  the  "hard  sphere"  value  of  Ou  „  which  gives  the 

0  """a 

predicted  binary  diffusion  coefficient  at  300  K.  The  hard  sphere  diameter  of 

the  hydrogen  molecule  (whluii  glvos  the  correct  viscosity  of  pure  Hj  at  300°)  is 
0 

about  2.69^''  A.  Them 

orchard  sphere)  -  2o  (hard  sphere)  -  Om  (hard  sphere)  (109) 
n  ^2 

In  surveying  this  list  one  notices  that,  in  general,  the  smallest  hydrogen  atom 
diameters  are  associated  with  a  priori  estimates,  while  the  largest  diameters 
are  associated  with  experimental  determinations,  however  indirect.  Thus,  dif¬ 
fusion  coefficients  and,  hence  Lewis-Semenov  numbe’-'s,  based  on  these  a  priori 
estimates  will  exceed  considerably  those  based  on  the  work  of  Harteck,  Amdur 
and  Wise,  Inasmuch  as  this  diffusion  coefficient  Is  extremely  important  in 
many  diverse  applications^,  additional  experimental  data  (preferably,  using 


e.q, ,  laminar  flame  theory,  theory  of  explosion  limits,  heat  transfer  from 
combustion  gases,  loss  of  hydrogen  atoms  from  planetary  atmosphc'ces,  etc. 
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different  techniques)  are  needed  in  order  to  verify  the  accuracy  of  the  avail¬ 
able  data*  and  to  extend  the  temperature  range. 


Lewls-Semenov  Number  for  Hydrogen  Atom  Diffusion  in  Diatomic  Hydrogen  and 
Combustion  Products 

For  an  equilibrium  dissociating  diatomic  gas*  the  composition  (degree 
of  dissociation)  will  vary  rapidly  with  changes  in  pressure  and  temperature 
level.  While  this  change  in  composition  does  not  in  itself  appreciably  affect 
the  atom-molecule  binary  diffusion  coefficient,  it  does  affect  the  frozen  ther¬ 
mal  diffuslvlty  3vj/(pCp  j)  of  the  mixture  with  the  result  that  the  Lewis-Semenov 
number  le^.  decreases  from  values  in  excess  of  unity  (when  the  degree  of  disso¬ 
ciation  is  small)  to  values  lecc  than  unity  (when  the  gas  is  almost  completely 

dissociated).  To  illustrate  this  behavior,  wr  have  computed  values  of  Le- 

t  * 

corresponding  to  Uie  tabulated  results  of  Relsfeld  for  equilibrium  dissociating 
hydrogen  in  the  temperature  range  from  1CXX1°K  to  40C)0°K**,  at  pressures  of  0.1, 
1,  10,  and  100  atmospheres  (see  Fig. 16),  The  curves  labeled  "P  =0"  and 
"P  a  00"  merely  indicate  complete  dissociation  and  the  absence  of  dissociation, 
respectively,  for  all  teinpeiatuijs  between  1000°K  and  /lOOO^K.  When  hydrogen 
atoms  CO  exist  with  heavinr  molecules  than  (as  they  do  in  equilibrium  com¬ 
bustion  product  gases),  the  corresponding  values  of  Le^  are  usually  much  higher 
than  those  calculnled  above  for  pure  hydrogen.  To  demonstrate  this,  consider 
a  hypothetical  mixture  of  water  vapor  with  trace  amount?  of  hydrogen  atoms  at 

T  =  3000°K,  p  -  1  atmosphere.  The  diffusion  coefficient  D„  _  will  he  esti- 

H-HjO 

mated  by  combining  the  "hard  sphere"  values  of  cr„  and  n  at  3000  K  in  accord 

n  HjU 

with  the  rule: 


t 


For  this  calculation,  we  used  Reisfeld's  values  oft  C  for  pure  H  and  Hj,  the 

P 

equilibrium  mole  fraction  p  for  the  mixture  [from  v/hich  was  estimated 

using  Eq. (3.2-41)  of  reference  5l],  and  t.hp  binary  diffusion  coefficient  D„  „ 

H“n 
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Then  will  be  assumed  to  be  qiven  byj 

H-fljU 


i  i 

=  2.628  »  10->  (T  M  [(H„  +  Vo>''<2»m'«H,0>J  ("H-Hao) 


-3 


(111) 


Tl’ie  viscosity  of  water  vapor  at  this  temperature  is  first  estimated  from  the 

o  .  « 

Lennard-Jones  parameters  a  =  2.824  A,  e/K  =  230.9*  6  =  2.333  given  In  reference 
51.  The  hard  sphere  value  of  (T„  *  corresponding  to  this  viscosity  is  about 

O  HjU 

2,37  A.  Similarly,  the  hard  sphere  value  of  (J„  can  be  obtained  from  the  Harteck* 

H 

Amdur  formula  (l07)  combined  with  the  value  of  (j„  (hard  sphere).  As  above, 

”3  0 

a„  is  obtained  from  the  viscosity  of  pure  hydrogen  gas  at  3000  K.  Using  the 
da  0  A  0 

potential  parameters  j  =  2.915  A,  t./K  =  38  K  for  pure  hydrogen*  one  finds 

0  O  ’  0 

CJ„  2.27  A.  The  corresponding  value  of  o„  is  2.63  A.  Thus,  (Tu  „  -  =  2.50  A 
Ha  «  H  n-HjU 

1  «  _ 5/.  ^1 _ _ _ _ _ 1 _ _  -A _ A.  _ _ —  - _ 


The  Prandtl  number  Prj^  for  steam  at  tempera - 


and,  hence,  D„  „  -  -  50cm* (sec) 

H-HaO  Q 

tuies  in  the  range  3000  K  is  not  accurately  known,  but  is  probably  somewhat 
greater  than  the  lower  temperature  (800°K)  value®®  of  1,01.  Combining  this 
estimate  with  the  calculated  values  of  ^  and  ^  gives  a  Lewis-Semenov 
number  of  about  3.4.  For  the  stoichiometric  combustion  of  hydrogen  and  oxygen 
at  typical  rocket  motor  pressures,  the  major  constituents  of  the  product  gases 
are  steam,  diatomic  hydrogen,  hydroxyl  radicals  and  hydrogen  atoms,  with  steam 
contributing  something  of  the  order  of  70  percent  of  the  number  of  particles 
per  cubic  centimeter.  Tnua,  the  calculation  outlined  above  suggests  that  the 
Lewis-Semenov  number  for  hydrogen  atom  diffusion  through  such  mixtures^  is  in 


We  have  selected  hydrogen  atom  diffusion  as  an  example  becausei 

a.  using  tabulated  populations  and  boats  of  formation  of  the  three  prin¬ 
cipal  radicals  H,0,OH,  the  hydrogen  atoms  account  for  most  of  the  free 
stream  chemical  energy,  usually  more  than  50  percent 

b.  of  the  three  principal  radicals,  hydrogen  atoms  will  probably  diffuse 
most  rapidly  through  the  gas  mixture'*® 

c.  if  the  chemistry  in  the  gas  phase  is  such  that  radicals  can  reach  the 
confining  walls  of  clean  metal  hydrogen  atoms  will  probably  have  the 
largest  recombination  probability’** 
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the  range  of  3.5  even  when  conservative  estimates  are  introduced  for  the  dif¬ 
fusion  coefficient  q.  This  further  implies  that  heat  transfer  rates  from 
combustion  gases  (containing  appreciable  amounts  of  hydrogen  atoms)  to  high  tem¬ 
perature  solids  may  differ  appreciably  from  tnose  calculated  under  the  assump¬ 
tion  Le^  =  1.  In  the  case  of  oxy-acetylene  flame  gases  (at  one  atmosphere 
pressure)  this  writer  estimated  about  a  ?.0  percent  effect  on  the  heat  flux  at 
the  highest  surface  temperature  reported  (2000^^113  01(106)  by  Gledt,  Cnhb,  and 
Russ*’.  For  this  reason,  heat  flux  data  at  surface  temperatures  approaching 
the  flame  temperature  would  be  extremely  interesting. 
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Vm  CQNCLUDIMS  REMARKS 


We  have  attempted  to  outline  the  way  In  which  exieting  heat  tranefe:: 
data  and  physico-chemical  data  can  be  brought  to  bear  on  the  problem  of  pre¬ 
dicting  heat  transfer  rates  in  chemically  reacting  syeteroa.  The  methods 
suggested  herein  are  equivalent  to  the  following  set  of  hypotheses! 

1.  The  total  boat  flux  can  be  regarded  as  being  the  sum  of  a 
conductive^  (thermal)  contribution  and  a  diffusion-chemical 
reaction  contribution  each  calculated  as  if*  the  location 
of  the  chemical  change  were  confined  to  the  Interface. 

2.  The  conductive  contribution  to  the  heat  flux  is  then  the 
result  of  the  solid  and  gas  not  being  in  equilibrium  with 
respect  to  thermal  and  kinetic  energy  and  is  estimated  from 
the  equation! 

~  St(Ra,Prj^^^)  G  {Ah^  +  rj^  Ah,^j  J  (ll2) 

3.  The  contribution  Ajj  due  to  diffusion  and  subsequent  chemical 
change  is  estimated  from  the  equation! 

=ZVl 

i  i 

This  contribution  is  then  the  result  of  the  solid  not  being 
in  chemical  equilibrium  with  the  gas  mixture*'^**'. 

4.  The  non-dimensional  heat  and  mass  transfer  coefficients 
(Stanton  numbers)  may  be  estimated  from  existing  non-reactive 
data  or  theory  provided  some  account  is  taken  of  the  effect 
of  variable  fluid  properties, 

^Tn  this  context,  the  use  of  conductive  is  not  intended  to  convey  the  absence 
of  fluid  motion  (convection) 

*If  this  is  not  the  case,  then  each  term  is  in  error  but  their  sum  is  approxi¬ 
mately  unchanged. 


If  these  statements  are  combined*  as  done  in  the  text*  then  the  heat  transfer 
rate  will  be  given  byi 


^  G  Ah' 


{l  t  (r^  -  1)^  *  I  (r„  ,  -  1)^^} 

I  X  ^0  ^  D,1  ^0  J 


ClU) 


where  the  recovery  factors  r_  ,  for  chemical  energy  may  be  related  to  the 

1 

Reynolds  analogy  factors  a(Pr|Re)  byi 


St(Re,Prjj^^) 

'*0,1  ~  St(Re,Prj^*  j)  ~  sfPr^*  j|Re) 


(115) 


For  laminar  boundary  layer  flow,  the  Reynolds  analogy  factors  s(Pr)  are  Indepen 
dent  of  the  local  Reynolds  number.  Exact,  values  of  8(Pr)  for  the  case  of  con¬ 
stant  property  laminar  boundary  layer  flow  over  a  flat  plate  are  shown  plotted 
In  Fig.  17^.  If  both  Pr^^  ^  and  the  Prjj  ^  are  not  too  small  compared  with  unity 
It  will  be  found  thati 


'’D.i  * 


I 


(116) 


where  Pr.  ^/Pr^,  .  Is  recognized  as  the  Lewis-Semenov  number  Lo-  ,  for  the  diffus 

AttUyl  tfl 

Ion  of  species  1.  For  turbulent  boundary  layer  flow,  eq,(ll6)  la  often  a  good 
approximation.  Actually,  In  this  case,  the  Reynolds  analogy  factor  should  have 
a  weak  Reynolds  number  dependence.  This  Is  shown  in  Fig. 13  where  the  Reynolds 
analogy  factor  is  plotted  against  the  Prandtl  number^  for  several  values  of  the 


^Constructed  from  tabular  values  given  by  Eckert*’  and  Merk’-’ 

*This  figure  should  not  be  relied  on  for  Prandtl  numbers  which  are  very  differ¬ 
ent  from  unity,  since  the  seml-emplrica 1  theory  upon  which  it  is  based  breaks 
down  in  tliujitf  exLi'ttines 
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FIG.  17  REYNOLDS  ANALOGY  FACTORS  FOR  THE  LAMINAR  BOUNDARY 
LAYER  FLOW  OVER  A  FLAT  PLATE 


Reynolds  number  based  on  tube  diameter  (for  fully  developed  pipe  liow). 

It  is  observed  that  when  the  recovery  factors  r^^,  r^  j  are  equal  to 
unltyi  the  bracketed  term  in  eq.(ll4)  reduces  to  unity  and  the  total  enthalpy 
difference  Ah°  across  the  boundary  layer  becomes  the  true  driving  force  for 
energy  transport.  More  generally,  the  driving  force  will  bet 

\  “'kin  +  I 'd,!  “’ch.m.l 
1 

when  St(Re,Pr.  ,)  is  used  as  the  energy  transport  coefficient. 

A,  f 

The  conditions  under  which  the  terms » 

Ah 

(r_  ,  -  1)— (118) 

Ah° 

are  expected  to  be  appreciable  have  been  discussed  in  Section  III,  where  empha¬ 
sis  has  been  given  to  the  case  of  local  thermochemica 1  equilibrium.  Numerical 
examples  for  the  case  of  disaociatlnq  hydrogen  over  a  range  of  pressures  and 
temperatures  show  that  large  values  of  are  attained  wheni 

1.  the  temperature  difference  across  the  boundary  layer  is  smallj 

2.  the  temperature  level  is  in  the  range  of  maximum  chemical  con¬ 
tribution  to  the  heat  capacity  fl .e. .  maximum  c  .  /c  „)- 

-  PiChem  p,eq 

The  estimation  of  non-reactlve  film  conductances  for  the  case  of  turbulent 
boundary  layer  development  in  rocket  motor  nozzles  is  discussed  in  Section  IV. 

If  local  thermochemica 1  equilibrium  Is  achieved  in  the  gas  phase,  then  the  en¬ 
thalpy  difference  Ah'^  can  be  obtained  using  the  methods  illustrated  in  Section  V 
for  the  case  of  oxy-hydrogen  combustion.  Large  entha Ipy/mixture-ratio  charts 
have  been  included  for  this  system  at  three  total  pressures*  10,  30  and  60 
atmospheres.  In  general,  chemical  non-equilibrium  effects  strongly  influence 
the  establishment  of  Ah"  and  the  Ah  ,  ,/ilh’.  A  semi-quantltstive  discussion 

Cn0ni«  1 


of  the  Influence  of  chemical  kinetic-aerodynamic  parameters  on  the  heat  trans¬ 
fer  is  given  in  Section  VI.  The  functions  introduced  therein  are  knowni 
however «  in  only  a  very  limited  number  of  circumstances. 

Lastly!  the  estimation  of  the  transport  properties  of  gas  mixtures  has 
been  discussed  (Section  VI)  with  emphasis  on  the  diffusion  coefficient  for 
molecular  fragments.  Detailed  results  are  given  for  the  case  of  equilibrium 
dissociating  hydrogen  together  with  a  discussion  of  the  current  gaps  in  our 
knowledge  of  these  diffusion  coefficients. 

An  extensive  and  hopefully  representative  list  of  references  is  Included 
to  assist  the  reader  in  retracing  developments  which  have  been  compressed  because 
of  space  limitations.  Much  of  what  is  included  is  tentative,  intuitive  and 
unfinished.  It  is  hoped  that  these  speculations  and  gaps  will  be  clarified  by 
subsequent  theoretical  and  experimental  Investigations  in  this  highly  important 
discipline. 
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APPENDIX  1 

THERMODYNAMIC  PROPERTIES  OF  THE  HYCftOGEN/OXYGEN  SYSTEM 


A  large  number  of  equilibrium  calculations  exist  In  the  literature  for 
the  hydrogen/oxygen  system.  In  setting  out  to  construct  a  convenient  enthalpy- 
mixture  ratio  chart  at  constant  pressure,  it  was  our  original  Intention  to  make 
as  much  of  the  existing  work  as  possible.  It  soon  became  clear,  however,  that 
the  sources  differed  w.th  one  another  in  one  or  more  of  the  following  respects: 

1.  mixture  ratios  considered 

2.  pressures  considered 

3.  temperatures  considered 

4.  choice  of  eiiLiidlpy  basis 

5.  choice  of  input  data 

6.  manner  of  presentation  of  computational  results. 

A  review  of  this  work  led  to  the  conclusion  that  extensive  supplementary  cal¬ 
culations,  interpolation  and  readjustment  of  available  results  would  be  necessary 
to  make  available  even  one  convenient  enthalpy-mixture  ratio  chart  fl.e..  at 
one  pressure). 

As  a  result  the  desired  computations  were  progranwied  for,  and  carried 
out  by,  an  IBM  650  digital  machine  using  the  tabulated  thermochemical  data  of 
reference  56.  The  entha  Ipy/mlxture-ratio  results  at  the  three  selected  l.mal 
pressures  (lO,  30,  60  atmospheres)  are  included  in  the  form  of  large  graphs 
(see  foldouts  1,  2,  3  ),  Species  compositiens  corresponding  to  these  conditions 
are  shown  (to  reduced  scale)  in  the  text  (see  Figs.  10.  11.  12).  in  Table  4» 
we  list  several  references  in  which  graphical  or  tabular  data  bearing  on  the 
hydrogen/oxygen  system  are  given.  Most  of  the  references  treat  pure  hydrogen 
(<^  =oo  ),  pure  oxycen  (0  =0),  or  stoichimetrl c  hydrogen/oxygen  mixtures  (0  =  l). 


<1)  _  equivalence  ratio  -  (na/0j.)/(H2/02)g^Qj^j^  .  The  relation  between  4', f 
and  other  common  tr.ixture  ratio  parameters  is  given  in  Tabic  5 
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Useful  exceptions  are  the  work  of  Hottelj  WiUiams  anrf  Sat-terfiwlrf  (^949)5*, 
Reichert  (l950)**i  Herbert  and  Zlebland  (1958)^’,  Gordon  and  McBride  (1959)^* 
and  Baker  (1959)^.  The  pressure  levels  and  mixture  ratios  considered  in  each 
of  these  investigations  are  given  in  columns  2  and  3  (pressure  in  atmospheres, 
mixture  ratio  parameter  =  equivalence  ratio).  It  will  be  noted  that,  according 
to  Reichert  (l950)®*,  large  enthalpy /mixture-ratio  charts  for  the  hydrogen/oxy¬ 
gen  system^  at  the  total  pressures  0.1,  1.0,  10,  100  kg/cm**  can  be  obtained 
from  the  Ministry  of  Supply,  London.  The  "reaction  enthalpy"  which  appears  in 
the  work  of  I.utz  (1946)*’,  Reichert  (1950)**  and  Herbert  and  Zlebland  (1958)*’ 
is  nothing  but  the  total  (sensible  +  chemical)  enthalpy  of  the  gas  mixture,  the 
only  difference  being  the  choice  of  base  value  on  the  absolute  enthalpy  scale. 
This  accounts  for  the  negative  values  occurlng  in  the  tabular  and  graphical 
data  of  these  authors. 


Twenty  charts  were  completed  for  the  carbon/hydrogen/oxygen  system  at  the  total 
pressures  0.1,  1.0,  10,  100  kg/cm*i  i .e. ,  four  charts  for  each  of  the  follow¬ 
ing  carbon-to-hydrogen  ratiost  0,  0.2,  0.4,  0.6,  0.8 


*  1  kg/cm* 


=  .9678a tm 
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Reference  63»  for  hydrogen  and  steam,  was  not  a\’ailable  to  the  author  at  the  time  this  compilation 
was  prepared 
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TABLE  5 


CORRESPONDENCE  OF  MIXTURE-RATIO  PARAMETERS  IN  THE  HYDROGEN/OXYGEN  SYSTEM^ 


Equivalence 

Mass  Fraction 

Mass  Fraction 

Oa/Ha 

Ratio 

of  Oxygen  in 

of  Hydrogen  in 

Mass  Ratio 

Mixture 

Mixture 

<!• 

f 

1-f 

0.000 

l.OOOCi 

0.00000 

oo 

0.100 

0.9876 

0.0124i'f 

79.37 

0.200 

0.9754 

0.02458 

39.68 

0.300 

0.9636 

0.03642 

26.46 

0.400 

0.9520 

0.04798 

19.84 

0.500 

0.9407 

0.05927 

15.87 

0.600 

0.9297 

0.07029 

13.23 

0.700 

0.9189 

0.08105 

11.34 

0.800 

0.9084 

0.09157 

9.921 

0.900 

0.8981 

0.1019 

8.818 

1.000 

0.8881 

0.1119 

7.937 

1-100 

0.8783 

0.1217 

7.215 

1.200 

0.8637 

0.1313 

6.614 

1.300 

0.8593 

0.1407 

6. 105 

1.400 

0.8501 

0.1499 

5.669 

1.500 

0.8410 

0.1590 

5.291 

1.600 

0.8322 

0,1678 

4.960 

1.700 

0.8236 

0.1764 

4.669 

1 .800 

0.8151 

0.1849 

4.409 

1.90C' 

0.8068 

0.1932 

4.177 

2.000 

0.7987 

0.2013 

3.968 

2.500 

0.7605 

0.2395 

3.175 

3.000 

0.7257 

0.2743. 

2 . 646 

3.500 

0.6940 

0.3060 

2.268 

4.000 

0.6649 

0.3351 

1.984 

4.500 

0.6382 

0.3618 

1.764 

5.000 

0.6135 

0.3865 

1.587 

5.500 

0.5907 

0.4093 

1.443 

6.000 

0.5695 

0.4305 

1.323 

6.500 

0.5498 

0.4502 

3..221 

7.000 

0.5314 

0.4687 

1.134 

7.500 

0.5141 

0.4859 

1.058 

8.000 

0.4980 

0.5020 

0.9921 

8.500 

0.4829 

0.5171 

0.9337 

9.000 

0.4686 

0.5314 

0.8818 

9.500 

0.4552 

0.5448 

0.3354 

^See  footnote  on 

following  page 
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TABLE  5  (COigiNUED) 

CORRESPONDENCE  OF  MIXTURE-RATIO  PARAMETERS  IN  THE  HYDROGEN/OXYGEN  SYSTb/ 


Equiva Isnce 
Ratio 

Mass  Fraction 
of  Oxyqen  in 
Mixture 

Mass  Fraction 
of  Hydrogen  in 
Mixture 

Oj/Ha 

Mass  Ratio 

<I> 

f 

1-f 

r 

10.000 

0.4425 

0.5575 

0.7937 

20.000 

0.2841 

0.7159 

0.3968 

30.000 

0.2092 

0.7908 

0.2646 

40.000 

0.1656 

0.83U 

0.1984 

50.000 

0.1370 

0.8630 

0.1587 

oo 

0.0000 

1.0000 

0.0000 

t 


Calculated  from  tha  relations t 


-1 

f  =  {l  -t-  [  (4.032) /(32.00)]  <!>} 


-  32.00 
‘  ~  4.032  <t 
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TABLE  6 

PRODUCT  GAS  COMPOSITION^  AND  MEAN  MOLECULAR  WEIGHT 
FOR  STOICHIOMETRIC  OXV-HYDROGEN  C0MBUSTIC!NtP=10  ATMOSPHERES 


T(°K) 

Mole 

Fraction,  x 

HjO 

Ha 

Oa 

OH 

0 

H 

"  ^ 

1000 

.9990 

.0,1000 

.0^,7602 

.0io3601 

.Ot,l66l 

,0^02269 

18.0 

1200 

.9990 

.0,1000 

.0,ol601 

.0,1192 

.0,,3717 

.081971 

18.0 

1400 

.9990 

.0, 1000 

.0,2022 

.067595 

.0,4856 

18.0 

1600 

.9989 

.0,9999 

.O44368 

.0^5451 

.061077 

.0,1716 

18.0 

IfiOO 

.9985 

.0,9995 

.0,2892 

.0,1962 

.067252 

.0,3571 

18.0 

2000 

.9952 

.0,3089 

.0,8682 

.0,7804 

.0,6832 

. 0^28 51 

18.0 

2200 

.9376 

.0,7022 

.0,2609 

.0,2534 

.O44744 

•0, 1488 

17.9 

2400 

.9721 

.01484 

.0,5669 

,0,6496 

.0,2226 

.0,6109 

17.8 

2600 

.94U 

.02789 

.01067 

.01421 

.0,8145 

.0,2021 

17.7 

2A00 

.8993 

.04834 

,01726 

.02705 

.0,2403 

.0,5669 

17.4 

3000 

.8323 

.07430 

.02652 

.04658 

.0,6181 

.0136] 

16.9 

3200 

.7410 

.1030 

.03623 

.07200 

.01369 

.02909 

16.2 

3400 

.6270 

.1423 

.04629 

.1016 

.02721 

.05557 

15.3 

3600 

.4976 

.1732 

.05388 

.1302 

.04850 

.09649 

14.1 

3800 

.3654 

.1936 

.05750 

.1518 

.07853 

.1531 

12.8 

4000 

.2455 

.1978 

.05621 

.1609 

.1164 

.2232 

11.4 

The 

notation 

.0, 1000 

impliec  the 

number  .OOlOOOj  similar 

ly  the  notation 

•  ^inpllGS  tho  niiTnbor 
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TABLE  7 


PRODUCT  GAS  COMPOSITION  AND  MEAN  MOLECULAR  WEIGHT 
FOR  STOICHIOMETRIC  OXY-HYDROGEN  C0MBUSTI0N«P=30  ATMOSPHERES 


T(°K) 

Mole 

Fraction,  x 

HjO 

0, 

OH 

0 

H 

1000 

.9990 

.0,1000 

.0i ,2534 

.0,o2079 

.0,,553a 

.O10I3IO 

18.0 

1200 

.9990 

.0,1000 

.0,^5336 

.0,6881 

.0,,1239 

.0,1138 

18.0 

1400 

.9990 

.0,1000 

.0,6741 

.0,4385 

.0,ol623 

.0.,2803 

18.0 

1600 

.9990 

.0,1000 

.0,1456 

.0,9953 

.0,3591 

.0,3132 

18.0 

1800 

.9988 

.0,9998 

. 0^9640 

.0,1133 

.0,2417 

.0,2062 

18.0 

2000 

.9967 

.0, 1866 

.0,7960 

.0,5806 

.0,3777 

.0,1279 

18.0 

2200 

.9914 

.0,4537 

.0,2054 

.0,1817 

.O42430 

.0,6946 

18.0 

2400 

.9808 

.0,9974 

.0,4257 

.0,4615 

.0,1113 

.0,2892 

17.9 

2600 

.9616 

.01985 

.0,7286 

.0,9906 

.0,3885 

.0,9842 

17.8 

2800 

.9304 

.03395 

.01248 

.01928 

.n,i 180 

.0,2743 

17.6 

3000 

.8836 

.05447 

.01879 

.03346 

.0,3004 

.0,6704 

17.2 

32.00 

.8186 

.07975 

.02701 

.05343 

.0,6824 

.01443 

16.8 

3400 

.7345 

.non 

.03548 

.07817 

.01375 

.02820 

16.2 

3600 

.6334 

.1409 

.04402 

.1061 

.02531 

.05023 

15.4 

3800 

.5209 

.1696 

.05070 

.1335 

.04257 

.08275 

14.4 

4000 

.4052 

.1909 

.05483 

.1561 

.06635 

.1266 

13.3 
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TABI,E  8 

PRODUCT  GAS  COMPOSITIOM  AND  MEAM  MOLECULAR  WEIGHT 
FOR  STOICHIOMETRIC  OXY»HYDROGEW  COMBUSTlONi  ATMOSPHERES 


T^K) 

Mole  Fraction,  : 

< 

HjO 

Ha 

Oa 

OH 

0 

H 

2000 

.9974 

.0,1621 

.0,5280 

.0,4408 

.0,2175 

.0,8431 

18.0 

2200 

.9932 

.0,3777 

.0,1520 

.0,1418 

.0^1478 

.0+4457 

18.0 

2400 

.9848 

.0,8095 

.0,3257 

.0,3637 

..0+6888 

.0,1842 

17.9 

2600 

.9696 

.01574 

.0,5891 

.0,7931 

.0,2470 

.0,6197 

17.8 

2800 

.9448 

.02753 

.0,9789 

.01537 

.0,7387 

.0,1747 

17.7 

3000 

.9076 

.04403 

.01517 

.02703 

.0,1908 

.0,4262 

17.4 

3200 

.8555 

.06554 

.02184 

.04355 

.0,4339 

.0,9252 

17.1 

3400 

.7871 

.09166 

.02930 

.06487 

.0,8837 

.01821 

16.6 

3600 

.7031 

.1196 

.03758 

.09039 

.01654 

.03274 

15.9 

3800 

.606? 

.  1488 

.07,464 

.1173 

.02825 

.05430 

15.2 

4000 

.5018 

.1746 

.05023 

.1429 

.04490 

.08562 

14.2 
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UNCLASSIFIED  I  UNCLASSIFIED 


